(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[    117608,       3199]
NotebookOptionsPosition[    103903,       2831]
NotebookOutlinePosition[    104300,       2848]
CellTagsIndexPosition[    104257,       2845]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", "\"\<@\>\"", "]"}], ";"}]], "Input",
 CellChangeTimes->{3.507197762859375*^9, 3.5075432984375*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Clear", "::", "wrsym"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Symbol \[NoBreak]\\!\\(\[FormalA]\\)\[NoBreak] is \
Protected. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\
\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/wrsym\\\", ButtonNote -> \
\\\"Clear::wrsym\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5085439238125*^9, 3.50864616446875*^9, 
  3.508714490359375*^9, 3.508800551109375*^9, 3.50882265946875*^9, 
  3.50888420228125*^9, 3.5091674949375*^9, 3.509596472640625*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Clear", "::", "wrsym"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Symbol \[NoBreak]\\!\\(\[FormalB]\\)\[NoBreak] is \
Protected. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\
\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/wrsym\\\", ButtonNote -> \
\\\"Clear::wrsym\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5085439238125*^9, 3.50864616446875*^9, 
  3.508714490359375*^9, 3.508800551109375*^9, 3.50882265946875*^9, 
  3.50888420228125*^9, 3.5091674949375*^9, 3.50959647265625*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Clear", "::", "wrsym"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Symbol \[NoBreak]\\!\\(\[FormalC]\\)\[NoBreak] is \
Protected. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\
\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/wrsym\\\", ButtonNote -> \
\\\"Clear::wrsym\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5085439238125*^9, 3.50864616446875*^9, 
  3.508714490359375*^9, 3.508800551109375*^9, 3.50882265946875*^9, 
  3.50888420228125*^9, 3.5091674949375*^9, 3.50959647265625*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Clear :: \
wrsym\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \
calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \
\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5085439238125*^9, 3.50864616446875*^9, 
  3.508714490359375*^9, 3.508800551109375*^9, 3.50882265946875*^9, 
  3.50888420228125*^9, 3.5091674949375*^9, 3.509596472703125*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"<<", "data`"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"<<", "Matcher`"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"<<", "classifier`"}], ";"}]}], "Input",
 CellChangeTimes->{{3.50650087121875*^9, 3.50650087740625*^9}, 
   3.506569846625*^9, {3.507543291265625*^9, 3.50754330975*^9}, {
   3.507543351125*^9, 3.507543353703125*^9}, {3.508543838328125*^9, 
   3.50854383865625*^9}, {3.508543930046875*^9, 3.508543932578125*^9}, {
   3.508579601546875*^9, 3.5085796056875*^9}, {3.5085818531875*^9, 
   3.508581860671875*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"so", " ", "first"}], ",", " ", 
   RowBox[{
   "we", " ", "just", " ", "use", " ", "the", " ", "track", " ", "of", " ", 
    "matrix"}], ",", " ", 
   RowBox[{"it", " ", "is", " ", "J1"}]}], "*)"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.50864776196875*^9, 3.50864779565625*^9}, 
   3.508651382609375*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"BeginPackage", "[", "\"\<W`PCA`\>\"", "]"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"J1", "::", "usage"}], "=", "\"\<J1\>\""}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"mPCA", "::", "usage"}], "=", "\"\<\>\""}], ";"}], 
  "\[IndentingNewLine]"}], "\n", 
 RowBox[{
  RowBox[{"Begin", "[", "\"\<`PP`\>\"", "]"}], ";"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 InitializationCell->True,
 CellChangeTimes->{{3.50864785875*^9, 3.50864789284375*^9}, {
   3.508648862703125*^9, 3.508648878046875*^9}, {3.50864891303125*^9, 
   3.508648960609375*^9}, {3.50864916846875*^9, 3.50864924625*^9}, {
   3.508649328296875*^9, 3.5086493678125*^9}, {3.508649869796875*^9, 
   3.50864987021875*^9}, {3.508649961125*^9, 3.508649991765625*^9}, {
   3.508650069578125*^9, 3.50865029209375*^9}, {3.508650345*^9, 
   3.508650346953125*^9}, {3.50865040203125*^9, 3.5086504021875*^9}, {
   3.50865044075*^9, 3.508650575453125*^9}, {3.508651178890625*^9, 
   3.508651205015625*^9}, 3.508651382609375*^9, {3.50890690559375*^9, 
   3.5089069411875*^9}, {3.508907013671875*^9, 3.508907233953125*^9}, {
   3.508907268453125*^9, 3.50890730671875*^9}, {3.508907684140625*^9, 
   3.50890769446875*^9}, {3.509167517546875*^9, 3.50916751940625*^9}, {
   3.509167923578125*^9, 3.50916797015625*^9}, {3.5091683784375*^9, 
   3.50916839584375*^9}, {3.509594335140625*^9, 3.509594358203125*^9}, {
   3.509595118328125*^9, 3.509595121421875*^9}, {3.5095951563125*^9, 
   3.509595159265625*^9}, {3.509595194265625*^9, 3.509595201*^9}, {
   3.5095956054375*^9, 3.50959567615625*^9}, {3.50959576740625*^9, 
   3.50959578025*^9}, {3.509595895125*^9, 3.509595927234375*^9}, 
   3.509595959359375*^9, {3.509595995453125*^9, 3.50959599940625*^9}, {
   3.5095967845625*^9, 3.509596842234375*^9}, {3.510031656078125*^9, 
   3.510031658*^9}, {3.5100317349375*^9, 3.510031750390625*^9}, {
   3.510041763375*^9, 3.5100417641875*^9}, {3.510389130296875*^9, 
   3.510389233109375*^9}, 3.51072306128125*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "ss", " ", "has", " ", "dim", " ", "of", " ", "m", "*", "n", "*", "l"}], 
    ",", " ", 
    RowBox[{
    "m", " ", "classes", " ", "and", " ", "n", " ", "samples", " ", "each"}], 
    ",", " ", 
    RowBox[{"each", " ", "sample", " ", "is", " ", "l", " ", "vector"}]}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"J1", "[", "ss_", "]"}], ":=", 
     RowBox[{"Module", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"Sw", ",", "Sb", ",", "m", ",", "mm"}], "}"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"mm", "=", 
         RowBox[{"Map", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Mean", "[", "#", "]"}], "&"}], ",", "ss"}], "]"}]}], ";",
         "\[IndentingNewLine]", 
        RowBox[{"Sw", "=", 
         RowBox[{"Total", "[", 
          RowBox[{"Map", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"Covariance", "[", "#", "]"}], "&"}], ",", "ss"}], "]"}],
           "]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Length", "[", "Sw", "]"}], ">", "0"}], ",", 
          RowBox[{"Sw", "=", 
           RowBox[{"Tr", "[", "Sw", "]"}]}]}], "]"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"m", "=", 
         RowBox[{"Mean", "[", 
          RowBox[{"Flatten", "[", 
           RowBox[{"ss", ",", "1"}], "]"}], "]"}]}], ";", 
        "\[IndentingNewLine]", "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Length", "[", "m", "]"}], ">", "0"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Sb", "=", 
           RowBox[{"Total", "[", 
            RowBox[{"Map", "[", 
             RowBox[{
              RowBox[{
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"#", "-", "m"}], ")"}], ".", 
                RowBox[{"(", 
                 RowBox[{"#", "-", "m"}], ")"}]}], "&"}], ",", "mm"}], "]"}], 
            "]"}]}], ",", "\[IndentingNewLine]", 
          RowBox[{"Sb", "=", 
           RowBox[{"Total", "[", 
            RowBox[{"Map", "[", 
             RowBox[{
              RowBox[{
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"#", "-", "m"}], ")"}], "*", 
                RowBox[{"(", 
                 RowBox[{"#", "-", "m"}], ")"}]}], "&"}], ",", "mm"}], "]"}], 
            "]"}]}]}], "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"Sb", "/", "Sw"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"outpou", ":", " ", "Eigenvectors"}], ",", " ", "Eigenvalues", 
     " ", ",", " ", "pcs"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"mPCA", "[", "ss_", "]"}], ":=", 
     RowBox[{"Module", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
        "lambda", ",", "ev", ",", "ppc", ",", "psc", ",", "tvp", ",", "stm"}],
         "}"}], ",", "\[IndentingNewLine]", 
       StyleBox[
        RowBox[{
         StyleBox["(*",
          FontColor->RGBColor[0.5, 0, 0.5]], 
         StyleBox[" ",
          FontColor->RGBColor[0.5, 0, 0.5]], 
         RowBox[{
          RowBox[{
           RowBox[{
            StyleBox["this",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox[" ",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox["is",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox[" ",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox["how",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox[" ",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox["PCA",
             FontColor->RGBColor[0, 0, 1]], 
            StyleBox[" ",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox["works",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox[" ",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox["in",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox[" ",
             FontColor->RGBColor[0.5, 0, 0.5]], 
            StyleBox["math",
             FontColor->RGBColor[0.5, 0, 0.5]]}], 
           StyleBox["...",
            FontColor->RGBColor[0.5, 0, 0.5]]}], 
          StyleBox[".",
           FontColor->RGBColor[0.5, 0, 0.5]]}], 
         StyleBox["*)",
          FontColor->RGBColor[0.5, 0, 0.5]]}],
        FontSize->16], "\[IndentingNewLine]", 
       RowBox[{"(*", 
        RowBox[{"tvp", ",", " ", 
         RowBox[{"each", " ", "row", " ", "is", " ", "a", " ", "vector"}]}], 
        "*)"}], "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"{", 
          RowBox[{"psc", ",", "tvp"}], "}"}], "=", 
         RowBox[{"Chop", "[", 
          RowBox[{"KarhunenLoeveDecomposition", "[", 
           RowBox[{
            RowBox[{"ss", "\[Transpose]"}], ",", 
            RowBox[{"\"\<Centered\>\"", "\[Rule]", "True"}]}], "]"}], "]"}]}],
         ";", "\[IndentingNewLine]", 
        RowBox[{"psc", "=", 
         RowBox[{"psc", "\[Transpose]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"(*", 
         RowBox[{
          RowBox[{
           RowBox[{"Standardize", "[", 
            RowBox[{"ss", ",", "Mean", ",", 
             RowBox[{"1", "&"}]}], "]"}], ".", 
           RowBox[{"tvp", "\[Transpose]"}]}], "\[Equal]", "psc"}], "*)"}], 
        "\[IndentingNewLine]", 
        RowBox[{"lambda", "=", 
         RowBox[{"Chop", "[", 
          RowBox[{"Eigenvalues", "[", 
           RowBox[{"Covariance", "[", "ss", "]"}], "]"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"(*", " ", 
         RowBox[{
          RowBox[{"result", " ", "is", " ", "also", " ", "raw"}], "-", 
          "wised"}], "*)"}], "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"tvp", ",", "lambda", ",", "psc"}], "}"}]}]}], 
      "\[IndentingNewLine]", "]"}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 InitializationCell->True,
 CellChangeTimes->{{3.50864785875*^9, 3.50864789284375*^9}, {
   3.508648862703125*^9, 3.508648878046875*^9}, {3.50864891303125*^9, 
   3.508648960609375*^9}, {3.50864916846875*^9, 3.50864924625*^9}, {
   3.508649328296875*^9, 3.5086493678125*^9}, {3.508649869796875*^9, 
   3.50864987021875*^9}, {3.508649961125*^9, 3.508649991765625*^9}, {
   3.508650069578125*^9, 3.50865029209375*^9}, {3.508650345*^9, 
   3.508650346953125*^9}, {3.50865040203125*^9, 3.5086504021875*^9}, {
   3.50865044075*^9, 3.508650575453125*^9}, {3.508651178890625*^9, 
   3.508651205015625*^9}, 3.508651382609375*^9, {3.50890690559375*^9, 
   3.5089069411875*^9}, {3.508907013671875*^9, 3.508907233953125*^9}, {
   3.508907268453125*^9, 3.50890730671875*^9}, {3.508907684140625*^9, 
   3.50890769446875*^9}, {3.509167517546875*^9, 3.50916751940625*^9}, {
   3.509167923578125*^9, 3.50916797015625*^9}, {3.5091683784375*^9, 
   3.50916839584375*^9}, {3.509594335140625*^9, 3.509594358203125*^9}, {
   3.509595118328125*^9, 3.509595121421875*^9}, {3.5095951563125*^9, 
   3.509595159265625*^9}, {3.509595194265625*^9, 3.509595201*^9}, {
   3.5095956054375*^9, 3.50959567615625*^9}, {3.50959576740625*^9, 
   3.50959578025*^9}, {3.509595895125*^9, 3.509595927234375*^9}, 
   3.509595959359375*^9, {3.509595995453125*^9, 3.50959599940625*^9}, {
   3.5095967845625*^9, 3.509596842234375*^9}, {3.510031656078125*^9, 
   3.510031658*^9}, {3.5100317349375*^9, 3.510031750390625*^9}, {
   3.510041763375*^9, 3.5100417641875*^9}, {3.510389130296875*^9, 
   3.510389172*^9}, {3.510389207*^9, 3.51038920734375*^9}, {
   3.510389257484375*^9, 3.510389264453125*^9}, 3.510723166828125*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"End", "[", "]"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"EndPackage", "[", "]"}], ";"}], "\[IndentingNewLine]"}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 InitializationCell->True,
 CellChangeTimes->{{3.50864785875*^9, 3.50864789284375*^9}, {
   3.508648862703125*^9, 3.508648878046875*^9}, {3.50864891303125*^9, 
   3.508648960609375*^9}, {3.50864916846875*^9, 3.50864924625*^9}, {
   3.508649328296875*^9, 3.5086493678125*^9}, {3.508649869796875*^9, 
   3.50864987021875*^9}, {3.508649961125*^9, 3.508649991765625*^9}, {
   3.508650069578125*^9, 3.50865029209375*^9}, {3.508650345*^9, 
   3.508650346953125*^9}, {3.50865040203125*^9, 3.5086504021875*^9}, {
   3.50865044075*^9, 3.508650575453125*^9}, {3.508651178890625*^9, 
   3.508651205015625*^9}, 3.508651382609375*^9, {3.50890690559375*^9, 
   3.5089069411875*^9}, {3.508907013671875*^9, 3.508907233953125*^9}, {
   3.508907268453125*^9, 3.50890730671875*^9}, {3.508907684140625*^9, 
   3.50890769446875*^9}, {3.509167517546875*^9, 3.50916751940625*^9}, {
   3.509167923578125*^9, 3.50916797015625*^9}, {3.5091683784375*^9, 
   3.50916839584375*^9}, {3.509594335140625*^9, 3.509594358203125*^9}, {
   3.509595118328125*^9, 3.509595121421875*^9}, {3.5095951563125*^9, 
   3.509595159265625*^9}, {3.509595194265625*^9, 3.509595201*^9}, {
   3.5095956054375*^9, 3.50959567615625*^9}, {3.50959576740625*^9, 
   3.50959578025*^9}, {3.509595895125*^9, 3.509595927234375*^9}, 
   3.509595959359375*^9, {3.509595995453125*^9, 3.50959599940625*^9}, {
   3.5095967845625*^9, 3.509596842234375*^9}, {3.510031656078125*^9, 
   3.510031658*^9}, {3.5100317349375*^9, 3.510031750390625*^9}, {
   3.510041763375*^9, 3.5100417641875*^9}, {3.510389130296875*^9, 
   3.510389172*^9}, {3.510389207*^9, 3.51038920734375*^9}, {
   3.510389257484375*^9, 3.510389264453125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{"nxrData", "[", 
    RowBox[{"-", "1"}], "]"}]}], ";"}], "\n", 
 RowBox[{"Dimensions", "[", "data", "]"}]}], "Input",
 CellChangeTimes->{{3.50650087121875*^9, 3.50650087740625*^9}, 
   3.506569846625*^9, {3.507543291265625*^9, 3.50754330975*^9}, {
   3.507543351125*^9, 3.507543353703125*^9}, {3.508543838328125*^9, 
   3.50854383865625*^9}, {3.508543930046875*^9, 3.508543940765625*^9}, 
   3.508576904296875*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"50", ",", "2", ",", "6", ",", "46"}], "}"}]], "Output",
 CellChangeTimes->{{3.50854393903125*^9, 3.508543950234375*^9}, 
   3.50864617196875*^9, 3.50871449453125*^9, 3.5088005851875*^9, 
   3.508822662953125*^9, 3.50888421278125*^9, 3.50916752434375*^9, 
   3.509186791171875*^9, 3.509596477046875*^9, 3.510031018046875*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"fft", "=", 
   RowBox[{"Flatten", "[", 
    RowBox[{"data", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"1", ",", "2"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"3", ",", "4"}], "}"}]}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{"Dimensions", "[", "fft", "]"}]}], "Input",
 CellChangeTimes->{{3.507543453484375*^9, 3.507543482453125*^9}, {
   3.507543535484375*^9, 3.50754354509375*^9}, {3.507543616140625*^9, 
   3.507543630296875*^9}, {3.508543971015625*^9, 3.50854399078125*^9}, 
   3.508576909078125*^9, {3.509167545171875*^9, 3.5091675745625*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"100", ",", "276"}], "}"}]], "Output",
 CellChangeTimes->{
  3.508543994359375*^9, 3.50856537603125*^9, 3.50864617940625*^9, 
   3.508714496984375*^9, 3.508800586875*^9, 3.50882266428125*^9, 
   3.508884214484375*^9, {3.509167525671875*^9, 3.509167577703125*^9}, 
   3.509596478515625*^9, 3.510031021109375*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"each", " ", 
    StyleBox["row",
     FontColor->RGBColor[0.5, 0, 0.5]], " ", "is", " ", "a", " ", "sample"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"sample", "=", 
     RowBox[{"Transpose", "[", "fft", "]"}]}], ";"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"sample", "=", "fft"}], ";"}]}]], "Input",
 CellChangeTimes->{{3.5085608461875*^9, 3.5085608689375*^9}, {
  3.508562194484375*^9, 3.508562222421875*^9}, {3.51003102321875*^9, 
  3.510031025015625*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"Needs", "[", "\"\<MultivariateStatistics`\>\"", "]"}], ";"}], 
  "*)"}]], "Input",
 CellChangeTimes->{3.506502541703125*^9, 3.50854391328125*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Chop", "[", "sample", "]"}], "==", "sample"}]], "Input",
 CellChangeTimes->{{3.50856587325*^9, 3.50856591109375*^9}}],

Cell[BoxData["True"], "Output",
 CellChangeTimes->{{3.50856589053125*^9, 3.508565912203125*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Chop", "[", 
   RowBox[{"Mean", "[", "sample", "]"}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.50857268571875*^9, 3.5085727044375*^9}, {
  3.5085769286875*^9, 3.50857693284375*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", 
   ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
   "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]], "Output",
 CellChangeTimes->{3.50857693021875*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"Covariance", "[", "sample", "]"}], ";", "\[IndentingNewLine]", 
   RowBox[{"Dimensions", "[", "%", "]"}]}], "*)"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.508560965796875*^9, 3.5085610034375*^9}, 
   3.508562263328125*^9, 3.508714512203125*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"od", "=", 
    RowBox[{"Ordering", "[", "cv", "]"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"i", "=", "1"}], ";", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"While", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"Total", "[", 
         RowBox[{"cv", "[", 
          RowBox[{"[", 
           RowBox[{"Ordering", "[", 
            RowBox[{"cv", ",", 
             RowBox[{"-", "i"}]}], "]"}], "]"}], "]"}], "]"}], "/", 
        RowBox[{"Total", "[", "cv", "]"}]}], "<", "0.8"}], ",", 
      RowBox[{"i", "++"}]}], "]"}], "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Total", "[", 
      RowBox[{"cv", "[", 
       RowBox[{"[", 
        RowBox[{"Ordering", "[", 
         RowBox[{"cv", ",", 
          RowBox[{"-", "5"}]}], "]"}], "]"}], "]"}], "]"}], "/", 
     RowBox[{"Total", "[", "cv", "]"}]}]}]}], "\[IndentingNewLine]", 
  "*)"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.5065030771875*^9, 3.506503100171875*^9}, {
   3.50650313859375*^9, 3.506503144765625*^9}, {3.506503205171875*^9, 
   3.50650320778125*^9}, 3.506511614125*^9, {3.50754358490625*^9, 
   3.507543589046875*^9}, 3.507543641578125*^9, {3.50754368415625*^9, 
   3.507543692828125*^9}, {3.507543794640625*^9, 3.507543836046875*^9}, {
   3.50856087278125*^9, 3.508560873734375*^9}, 3.508561571109375*^9, {
   3.508561746875*^9, 3.508561760359375*^9}, {3.508562280109375*^9, 
   3.50856228175*^9}, {3.508562630515625*^9, 3.50856263178125*^9}, {
   3.508563183265625*^9, 3.508563283296875*^9}, {3.508564809265625*^9, 
   3.508564811328125*^9}, 3.50856489021875*^9, {3.50856503253125*^9, 
   3.5085650733125*^9}, {3.508565923125*^9, 3.5085659316875*^9}, {
   3.508573807484375*^9, 3.50857381115625*^9}, 3.508575729078125*^9, {
   3.508576274234375*^9, 3.508576335125*^9}, {3.508576680015625*^9, 
   3.508576697625*^9}, 3.508576806171875*^9, 3.5085769419375*^9, {
   3.5085774904375*^9, 3.508577492*^9}, 3.508714512203125*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"pc", "=", 
   RowBox[{"PrincipalComponents", " ", "[", "sample", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Dimensions", "[", "pc", "]"}]}], "Input",
 CellChangeTimes->{{3.5075437084375*^9, 3.507543728109375*^9}, {
   3.507543817640625*^9, 3.507543822765625*^9}, 3.507543866421875*^9, {
   3.5075440229375*^9, 3.507544025765625*^9}, {3.507544100140625*^9, 
   3.507544283078125*^9}, {3.507544398078125*^9, 3.507544403875*^9}, {
   3.5085596245*^9, 3.508559637796875*^9}, {3.50856149375*^9, 
   3.508561504359375*^9}, {3.508561563359375*^9, 3.508561586390625*^9}, {
   3.50856170440625*^9, 3.508561705453125*^9}, {3.508562637078125*^9, 
   3.508562644359375*^9}, {3.50857572*^9, 3.50857574415625*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"100", ",", "276"}], "}"}]], "Output",
 CellChangeTimes->{3.5087145299375*^9, 3.508800592765625*^9, 
  3.508822669359375*^9, 3.5088842193125*^9, 3.510031048203125*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"tm", ",", "l", ",", "p"}], "}"}], "=", 
   RowBox[{"mPCA", "[", "sample", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Chop", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Total", "[", 
       RowBox[{"pc", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "i"}], "]"}], "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"Length", "[", 
         RowBox[{"pc", "[", 
          RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "}"}]}], "]"}], "]"}], 
   "*)"}]}]}], "Input",
 CellChangeTimes->{{3.508578315359375*^9, 3.50857839175*^9}, {
  3.510031077484375*^9, 3.51003112278125*^9}, {3.510031191328125*^9, 
  3.5100312161875*^9}, {3.5100312619375*^9, 3.510031283578125*^9}, {
  3.51003134159375*^9, 3.51003136834375*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"only", " ", "99", " ", "nonzero", " ", "pcs"}], "*)"}]], "Input",
 CellChangeTimes->{{3.510031888375*^9, 3.510031896546875*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Length", "[", 
  RowBox[{"Select", "[", 
   RowBox[{"l", ",", 
    RowBox[{
     RowBox[{"#", "\[NotEqual]", "0"}], "&"}]}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.510031372203125*^9, 3.510031501484375*^9}, {
  3.510031537*^9, 3.5100315581875*^9}}],

Cell[BoxData["99"], "Output",
 CellChangeTimes->{{3.510031378078125*^9, 3.51003144353125*^9}, {
   3.510031480578125*^9, 3.510031497375*^9}, 3.510031559234375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"p", "[", 
  RowBox[{"[", 
   RowBox[{";;", ",", 
    RowBox[{"100", ";;", "102"}]}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.510031570328125*^9, 3.510031608296875*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.510031584375*^9, 3.510031609078125*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"tpc", "=", 
   RowBox[{"pc", "[", 
    RowBox[{"[", 
     RowBox[{"All", ",", 
      RowBox[{";;", "99"}]}], "]"}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.508578006359375*^9, 3.508578031640625*^9}, {
  3.50857840578125*^9, 3.5085784111875*^9}, {3.508579329734375*^9, 
  3.50857933053125*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"pps", "=", 
    RowBox[{"tpc", "[", 
     RowBox[{"[", 
      RowBox[{";;", ",", 
       RowBox[{";;", "3"}]}], "]"}], "]"}]}], ";"}], "*)"}]], "Input",
 CellChangeTimes->{{3.508581333359375*^9, 3.50858137603125*^9}, 
   3.508714554546875*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"ListPlot3D", "[", "pps", "]"}], "*)"}]], "Input",
 CellChangeTimes->{{3.50858138309375*^9, 3.50858142409375*^9}, {
   3.508581526984375*^9, 3.5085815310625*^9}, 3.508581697390625*^9, 
   3.508581736625*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.508581731609375*^9, 3.508581731640625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"pick", " ", "only", " ", "one", " ", "factor"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"p1", "=", 
     RowBox[{"tpc", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", "1"}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"sp1", "=", 
     RowBox[{"Partition", "[", 
      RowBox[{"p1", ",", "2"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"Dimensions", "[", "%", "]"}]}]}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.50857931578125*^9, 3.50857933996875*^9}, {
   3.508579401578125*^9, 3.50857940375*^9}, {3.5085800598125*^9, 
   3.508580090921875*^9}, 3.508714567375*^9, 3.50888379915625*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"50", ",", "2"}], "}"}]], "Output",
 CellChangeTimes->{3.508884224765625*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"ddiff", ",", "dsame"}], "}"}], "=", 
   RowBox[{"eucDist", "[", "sp1", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"ns", ",", "nd"}], "}"}], "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Length", "[", "dsame", "]"}], ",", 
     RowBox[{"Length", "[", "ddiff", "]"}]}], "}"}]}], ";"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.50857968046875*^9, 3.508579750046875*^9}, 
   3.50858011725*^9, 3.508714567375*^9}],

Cell[BoxData[{
 RowBox[{"Min", "[", "dsame", "]"}], "\[IndentingNewLine]", 
 RowBox[{"Max", "[", "dsame", "]"}], "\[IndentingNewLine]", 
 RowBox[{"mdsame", "=", 
  RowBox[{"Mean", "[", "dsame", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ddsame", "=", 
   RowBox[{"StandardDeviation", "[", "dsame", "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Min", "[", "ddiff", "]"}], "\[IndentingNewLine]", 
 RowBox[{"Max", "[", "ddiff", "]"}], "\[IndentingNewLine]", 
 RowBox[{"mddiff", "=", 
  RowBox[{"Mean", "[", "ddiff", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"dddiff", "=", 
  RowBox[{"StandardDeviation", "[", "ddiff", "]"}]}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{3.508714567375*^9}],

Cell[BoxData["0.05302526212729708`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{
  3.508579816796875*^9, {3.508714567375*^9, 3.508714591765625*^9}, 
   3.50871642190625*^9, 3.508822791671875*^9, 3.508884230125*^9}],

Cell[BoxData["19.28253920672684`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{
  3.508579816796875*^9, {3.508714567375*^9, 3.508714591765625*^9}, 
   3.50871642190625*^9, 3.508822791671875*^9, 3.508884230125*^9}],

Cell[BoxData["2.316462344143316`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{
  3.508579816796875*^9, {3.508714567375*^9, 3.508714591765625*^9}, 
   3.50871642190625*^9, 3.508822791671875*^9, 3.508884230140625*^9}],

Cell[BoxData["3.7384425681952047`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{
  3.508579816796875*^9, {3.508714567375*^9, 3.508714591765625*^9}, 
   3.50871642190625*^9, 3.508822791671875*^9, 3.508884230203125*^9}],

Cell[BoxData["0.00038705485321344213`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{
  3.508579816796875*^9, {3.508714567375*^9, 3.508714591765625*^9}, 
   3.50871642190625*^9, 3.508822791671875*^9, 3.508884230203125*^9}],

Cell[BoxData["72.47845698955447`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{
  3.508579816796875*^9, {3.508714567375*^9, 3.508714591765625*^9}, 
   3.50871642190625*^9, 3.508822791671875*^9, 3.508884230203125*^9}],

Cell[BoxData["16.33636341883875`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{
  3.508579816796875*^9, {3.508714567375*^9, 3.508714591765625*^9}, 
   3.50871642190625*^9, 3.508822791671875*^9, 3.508884230203125*^9}],

Cell[BoxData["13.338353099618269`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{
  3.508579816796875*^9, {3.508714567375*^9, 3.508714591765625*^9}, 
   3.50871642190625*^9, 3.508822791671875*^9, 3.508884230203125*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"dx", "=", "0.4"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dds", "=", 
   RowBox[{
    RowBox[{"BinCounts", "[", 
     RowBox[{"dsame", ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "60", ",", "dx"}], "}"}]}], "]"}], "/", "ns"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ddd", "=", 
    RowBox[{
     RowBox[{"BinCounts", "[", 
      RowBox[{"ddiff", ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "60", ",", "dx"}], "}"}]}], "]"}], "/", "nd"}]}], 
   ";"}], "\[IndentingNewLine]"}], "\n", 
 RowBox[{
  RowBox[{"dist", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"i", "*", "dx"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "0", ",", 
       RowBox[{
        RowBox[{"Length", "[", "ddd", "]"}], "-", "1"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pdds", "=", 
   RowBox[{"Transpose", "[", 
    RowBox[{"{", 
     RowBox[{"dist", ",", "dds"}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pddd", "=", 
   RowBox[{"Transpose", "[", 
    RowBox[{"{", 
     RowBox[{"dist", ",", "ddd"}], "}"}], "]"}]}], ";"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.50857984678125*^9, 3.508579848203125*^9}, {
   3.50857992215625*^9, 3.508579952953125*^9}, 3.5085800103125*^9, 
   3.508714591765625*^9}]
}, {11}]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"pdata", "=", 
   RowBox[{"sp1", "[", 
    RowBox[{"[", 
     RowBox[{"All", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"templ", "=", 
   RowBox[{"sp1", "[", 
    RowBox[{"[", 
     RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tg", "=", 
   RowBox[{"ConstantArray", "[", 
    RowBox[{
     RowBox[{"ConstantArray", "[", 
      RowBox[{"0", ",", 
       RowBox[{"Length", "[", "templ", "]"}]}], "]"}], ",", 
     RowBox[{"Length", "[", "pdata", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"i", "=", "1"}], ";", 
     RowBox[{"k", "=", "1"}], ";", 
     RowBox[{"n", "=", "1"}]}], ",", 
    RowBox[{"i", "\[LessEqual]", 
     RowBox[{"Length", "[", "pdata", "]"}]}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"m", "=", "1"}], ",", 
       RowBox[{"m", "\[LessEqual]", "1"}], ",", 
       RowBox[{"m", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"tg", "[", 
          RowBox[{"[", 
           RowBox[{"n", ",", "k"}], "]"}], "]"}], "=", "1"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"n", "++"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"k", "++"}], ";"}]}], "\[IndentingNewLine]", "]"}], 
  ";"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{3.508883855265625*^9}],

Cell[BoxData[
 RowBox[{"EER", "[", 
  RowBox[{"pdata", ",", "templ", ",", "tg", ",", "classifyDist1", ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "10", ",", "0.01"}], "}"}]}], "]"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.50882296896875*^9, 3.508822997984375*^9}, {
   3.5088231206875*^9, 3.508823131859375*^9}, {3.508823209015625*^9, 
   3.50882321015625*^9}, {3.50883407290625*^9, 3.50883407815625*^9}, 
   3.508883855265625*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"3.289999999999974`", ",", 
   RowBox[{"{", 
    RowBox[{"0.14816326530612245`", ",", "0.14`"}], "}"}]}], "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{
  3.508823028609375*^9, 3.508823153390625*^9, 3.508823251109375*^9, 
   3.508833759609375*^9, 3.508833869625*^9, 3.508834029296875*^9, 
   3.508834100578125*^9, {3.50883417709375*^9, 3.508834208515625*^9}, 
   3.508883855265625*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"pg", "=", 
   RowBox[{"classifyDist1", "[", 
    RowBox[{"pdata", ",", "templ", ",", "3.282"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"VF", "[", 
  RowBox[{"tg", ",", "pg"}], "]"}]}], "Input",
 CellChangeTimes->{{3.508822878921875*^9, 3.50882287934375*^9}, 
   3.508822958609375*^9, 3.5088232355625*^9, {3.508824555859375*^9, 
   3.50882459328125*^9}, {3.508824798125*^9, 3.50882482084375*^9}, {
   3.5088333760625*^9, 3.5088334009375*^9}, {3.50883412171875*^9, 
   3.508834123859375*^9}, {3.508834206859375*^9, 3.508834212125*^9}, 
   3.508883943625*^9, {3.508884247890625*^9, 3.50888426265625*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.1484`", ",", "0.14857142857142858`", ",", "0.14`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.50882480059375*^9, 3.50883181859375*^9, 3.5088324579375*^9, 
   3.508832793921875*^9, {3.508833376703125*^9, 3.508833402*^9}, {
   3.508834119765625*^9, 3.50883413409375*^9}, {3.508884242640625*^9, 
   3.508884263015625*^9}}]
}, Closed]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"mf", "[", "x_", "]"}], ":=", 
   RowBox[{
    RowBox[{"VF", "[", 
     RowBox[{"tg", ",", 
      RowBox[{"classifyDist1", "[", 
       RowBox[{"pdata", ",", "templ", ",", "x"}], "]"}]}], "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.50883463696875*^9, 3.5088346904375*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"mf", "[", "1.5", "]"}]], "Input",
 CellChangeTimes->{{3.508834724109375*^9, 3.50883473375*^9}}],

Cell[BoxData["0.074`"], "Output",
 CellChangeTimes->{{3.50883473075*^9, 3.50883473409375*^9}, 
   3.50888427190625*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
     RowBox[{
     "as", " ", "there", " ", "are", " ", "few", " ", "samples", " ", "in", 
      " ", "one", " ", "class"}], ",", 
     RowBox[{
     "it", " ", "is", " ", "hard", " ", "to", " ", "use", " ", "error", " ", 
      "rate", " ", "here"}]}],
    FontColor->RGBColor[0.5, 0, 0.5]], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Minimize", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", 
        RowBox[{"mf", "[", "x", "]"}]}], ",", 
       RowBox[{
        RowBox[{"x", ">", "1"}], "&&", 
        RowBox[{"x", "\[Element]", "Reals"}]}]}], "}"}], ",", "x"}], "]"}], 
   "*)"}]}]], "Input",
 CellChangeTimes->{{3.50882475109375*^9, 3.50882477709375*^9}, {
  3.50882490746875*^9, 3.508824937828125*^9}, {3.508834269703125*^9, 
  3.5088343989375*^9}, {3.5088344674375*^9, 3.50883449375*^9}, {
  3.508834581671875*^9, 3.5088346215625*^9}, {3.508834697421875*^9, 
  3.508834781546875*^9}, {3.508834823171875*^9, 3.50883483584375*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "0.02`"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", "\[Rule]", "1.`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.508822889265625*^9, 3.50882325115625*^9, {3.508824556875*^9, 
   3.508824593625*^9}, 3.508824788515625*^9, {3.50883425534375*^9, 
   3.5088342773125*^9}, {3.508834314765625*^9, 3.5088343993125*^9}, 
   3.5088344695625*^9, {3.5088346001875*^9, 3.5088346225625*^9}, {
   3.508834700140625*^9, 3.508834769828125*^9}, {3.508834830171875*^9, 
   3.5088348328125*^9}}]
}, Closed]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"feer", "[", 
      RowBox[{"pdata_", ",", "templ_", ",", "tg_"}], "]"}], ":=", 
     RowBox[{"Module", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"rr", ",", "de", ",", 
         RowBox[{"tde", "=", "1"}], ",", "i", ",", "mr"}], "}"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"For", "[", 
         RowBox[{
          RowBox[{"i", "=", "0"}], ",", 
          RowBox[{"i", "<", "5"}], ",", 
          RowBox[{"i", "=", 
           RowBox[{"i", "+", "0.005"}]}], ",", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"rr", "=", 
            RowBox[{
             RowBox[{"VF", "[", 
              RowBox[{"tg", ",", 
               RowBox[{"classifyDist1", "[", 
                RowBox[{"pdata", ",", "templ", ",", "i"}], "]"}]}], "]"}], 
             "[", 
             RowBox[{"[", 
              RowBox[{"{", 
               RowBox[{"2", ",", "3"}], "}"}], "]"}], "]"}]}], ";", 
           "\[IndentingNewLine]", 
           RowBox[{"de", "=", 
            RowBox[{"Abs", "[", 
             RowBox[{
              RowBox[{"rr", "[", 
               RowBox[{"[", "1", "]"}], "]"}], "-", 
              RowBox[{"rr", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", 
           "\[IndentingNewLine]", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"de", "<=", "tde"}], ",", 
             RowBox[{
              RowBox[{"tde", "=", "de"}], ";", 
              RowBox[{"mr", "=", "rr"}]}], ",", 
             RowBox[{"Break", "[", "]"}]}], "]"}]}]}], "]"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mr", ",", "tde"}], "}"}]}]}], "\[IndentingNewLine]", 
      "]"}]}], ";"}], "*)"}], "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.50882469390625*^9, 3.5088246954375*^9}, 
   3.50882490334375*^9, {3.50883161815625*^9, 3.508831794921875*^9}, {
   3.50883247203125*^9, 3.508832562109375*^9}, {3.508832595265625*^9, 
   3.508832605453125*^9}, {3.5088327715625*^9, 3.508832841390625*^9}, {
   3.508832992296875*^9, 3.508833289125*^9}, {3.5088333448125*^9, 
   3.50883334709375*^9}, {3.50883340984375*^9, 3.508833414859375*^9}, {
   3.508833464234375*^9, 3.50883352571875*^9}, 3.508834220734375*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"pdds", ",", "pddd"}], "}"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "60"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0.3"}], "}"}]}], "}"}]}]}], "]"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.50857988015625*^9, 3.508579907421875*^9}, {
   3.508579964109375*^9, 3.508579992390625*^9}, 3.508714598421875*^9}],

Cell[BoxData[
 GraphicsBox[{
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJxd0r1LVXEAh/GLQiA4GEkohWhvvmVZ195fzlXLl9JB56YgHdtrdnaJoOZm
pwZbGprtD3B0FEdpaNPjc/w98P3C5fJwzvn8Dpc78v7j+oeuVqv19fRTf59v
d2jraG33d/X9W739avzd9p9/Y3tNH5aeadc7rj737ywefPrZXP/ffP+qGq6D
t9Pc393h/i/V5ka9C6V5rqfpH9XfvXq95TpcX+f8Pbn/Ymn8S6V5vr80510u
zQbi/MHwroR3Nbyh8IaLx/3D8b4j8fy1cj/nXQ/vRpxv490M71a8/2h4Y+HZ
eOPhTYQ3Gd7t8Gy8qfDuhHc3vOm4fzp+T6+fPd6+F+fZ7H54Np7N+7Tj/9AO
byY8G8/GexDvZ7OH4dl4Nt6j8Gz2ODwbz8Z7Ep7NnoZn49l4z8Kz2fPwbDwb
70V4NnsZno1n41Xh2axTuvk7hG/zvM15NpsNz8az8Ww8m82FZ+PZeDaezeZL
8/vYeDaejWezV+HZeDaejWez1+HZeDaejWezhfBsPBvPxrPZYng2no1n49ls
KTwbz8az8Wy2HJ6NZ+PZeDZ7E56NZ+PZeDZ7G56NZ+PZeDZbCc/Gs/FsPJut
hmfj2Xg2nn0CjPfsgQ==
     "]]}, 
   {Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJxd1nlI1EEUB3ApEoKgwiilA03NNPNY19v8vb11d42KoBIq8QyT7LKy+7aw
S82iogMSo5v+iMiKDisKKkpKs8Qs84jynyK6IHJ9399bmmEXGWbeZ957M4pB
uaWzCgb5+Pgc6v96fupjy8xzjzufVWtHj3hGq/boeOiX7msHMP+knWj4/Dq3
tEYzxnnGN5nz+i+t7kdPdkFjlQaO1o265Hi/9iD2D6a/hgu7PtTVaEWFnuFL
Qy/W7elx6/5QChvdfayw74D29IlnDJN15kbQws47t14278H+kZTXu6ghP6oa
vh/NP1+SlPumBvGjKD/wltYSUYnzRlPwyPSTzbWV8Pzp//oCqMu36E/B3Cp4
Y+mc8WNb0Rg9n3HUXPHw9/O5++FNoKDWTwH5R6rhBUq+vD+Qatc4sp837ocf
RPPq76c2bd2H+Il0Lye8+GOenn+w+OyFSL7shdCVs56EKuCFUvH863dbN+6G
N4kGyvfbDS+Mlh69amlrqoA3mcparjxon7QT3mTFC6fkIdFnOrq2w4ugn98v
+3bM2AFvCtk2hZ3v2rsFXqTEsxcpPntTZT97UeKzF614MbiPzdgfo/gxdK07
pvDdjR3aQHhcrOTH58WKz54B70+vxyDx7BnkPM4nDu9zG7w48j9sO/aibyM8
I9VHp3W03F4Lz4j8VsMzUlRlwum309bAi8d72AwvngauY4zuJdCpr9qCV4vX
w0ug8Okp/U+qHF4C6lkJL1HOYy9RfPaSUF8ZvCRZZy9JzmMvmYyPDKHtVavg
JSv5pdDOffWu7on6/aXIfbHnXWcvlTTzlPDO5BXwUvF+18FLo/EvTf2fDfDS
JJ69NLzfcnjTlPy8c/bS5Tz20ulm8NTe98NXwktX7kOT/rCnoR8l8Aj3W4z3
Rfj9zIVPeC9LEU8Sz+eR4pnQzyXwTPCWwzPJfvZM6O8yeCb8/SiFZ5Z+smeW
ePbMEs+eWfJnz4z6C+BZpF7uj0Wp14L3VgTPgv7Ohqevz4Fnlf3sWWWdPSv6
sQCeVfJhz4r+zoBnk3j2bFjPgWeTfNmzST3s2ZBvPjw71rPh2ZV67cp92CVf
9uxSP3sOJT+HUq93zp5D6mPPgXyc8DJkzp4+T4GnzyPgefezl6F4mRLPXqac
z553nb1MxfOus+dE/VZ4TsVzSr/Zcyr1euPZc8k6ey7Fc0k+7LmU/FyK51bq
dUu/2HMr/XOT/n8Pe27Fy1K8LNnPXpbiZUk8e979/wD3dgMa
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  ImageSize->{552.5000000000002, Automatic},
  PlotRange->{{0, 60}, {0, 0.3}},
  PlotRangeClipping->True]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.508579889625*^9, 3.508579924890625*^9}, {
   3.508579955625*^9, 3.508580012875*^9}, {3.508714598421875*^9, 
   3.508714604609375*^9}, 3.508716512703125*^9, 3.50882279678125*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"use", " ", "a", " ", "few", " ", "more", " ", "factors"}], 
   "*)"}]}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.508580178015625*^9, 3.508580198546875*^9}, 
   3.508714647484375*^9, 3.508884079390625*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"ps", "=", 
   RowBox[{"tpc", "[", 
    RowBox[{"[", 
     RowBox[{";;", ",", 
      RowBox[{"1", ";;", "10"}]}], "]"}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"sps", "=", 
   RowBox[{"Partition", "[", 
    RowBox[{"ps", ",", "2"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Dimensions", "[", "%", "]"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.508580204921875*^9, 3.50858022*^9}, 
   3.508714647484375*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"50", ",", "2", ",", "10"}], "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{3.508580221578125*^9, 3.508714647484375*^9, 
  3.508884296140625*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"ddiff", ",", "dsame"}], "}"}], "=", 
   RowBox[{"eucDist", "[", "sps", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"ns", ",", "nd"}], "}"}], "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Length", "[", "dsame", "]"}], ",", 
     RowBox[{"Length", "[", "ddiff", "]"}]}], "}"}]}], ";"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.508580835375*^9, 3.508580837*^9}, 3.508714656546875*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dx", "=", "0.4"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dds", "=", 
   RowBox[{
    RowBox[{"BinCounts", "[", 
     RowBox[{"dsame", ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "60", ",", "dx"}], "}"}]}], "]"}], "/", "ns"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ddd", "=", 
    RowBox[{
     RowBox[{"BinCounts", "[", 
      RowBox[{"ddiff", ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "60", ",", "dx"}], "}"}]}], "]"}], "/", "nd"}]}], 
   ";"}], "\[IndentingNewLine]"}], "\n", 
 RowBox[{
  RowBox[{"dist", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"i", "*", "dx"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "0", ",", 
       RowBox[{
        RowBox[{"Length", "[", "ddd", "]"}], "-", "1"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pdds", "=", 
   RowBox[{"Transpose", "[", 
    RowBox[{"{", 
     RowBox[{"dist", ",", "dds"}], "}"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pddd", "=", 
   RowBox[{"Transpose", "[", 
    RowBox[{"{", 
     RowBox[{"dist", ",", "ddd"}], "}"}], "]"}]}], ";"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10002.},
 CellChangeTimes->{{3.508714656546875*^9, 3.50871469559375*^9}, 
   3.50871493403125*^9}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"pdds", ",", "pddd"}], "}"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "60"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0.2"}], "}"}]}], "}"}]}]}], "]"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10002.},
 CellChangeTimes->{
  3.50858086690625*^9, {3.508714656546875*^9, 3.50871469559375*^9}, 
   3.50871493403125*^9}],

Cell[BoxData[
 GraphicsBox[{
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJxd0rtqFVEARuGDQkCwUAyiKCHRROPd5CQn0UT3eI3XIj6CoJb2sfYNRDC1
dSqL2AgBS30AS0vRzspOx7XdC/8fwmFlZn8zCWfq8fNHT3YNBoPXf376z9zm
m35fyvbEy+/r21uF/lZ+rPe/2SkLw34/y4vxrbWvGx/r9V/1832pTPevuX93
x/3vyrOn/ca6/5+3p+N5H8rnT/321vvfVm9f3L+/Nf6B1pwfr+df1ecdbNfZ
oXYd73BrvCPhHW2NNxHeZDvP/ZPxvlNx/li8//HwpuP5Nt5MeCfi/U+GNxve
bPz9p8I7He93Jryz4dl458I7H96F8C7G/Ta+/ff4cC6eZ7P58Gy8+fh/DcMb
tuv1rcNbCM/GWwxvMbxReKPwRuEthWez5fi+L4dn410Kz2aX4/1sPBtvJTyb
rYZn49l4V8Kz2dXwbDwbr4Rns651/TqEb3Pe5nk2uxaejWfj2Xg2ux6ejWfj
2Xg2u9Ga/4+NZ+PZeDa7GZ6NZ+PZeDa7FZ6NZ+PZeDa7HZ6NZ+PZeDZbC8/G
s/FsPJvdCc/Gs/FsPJvdDc/Gs/FsPJvdC8/Gs/FsPJvdD8/Gs/FsPJs9CM/G
s/FsPJs9DM/Gs/FsPPs3AJj43Q==
     "]]}, 
   {Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJxd1mlIVFEUB3ApEoygRMmihdwyy1wnl3F5Z5ytWQyLwBDKsFSyUCEr08ol
NcUlDSwyoqAos40+lERFi30oqGgPS8owLSG/FFHRlxzPf87D+5hBDvfc3z3n
vePT4PzSdQXTfHx8jk58PT/Vq/u45xrQpsZjEhsSPNdP7Xboym+fZ1s0Xv+r
PdgcWfxlS743j3h9q8b50yllRsz5oZFyrajQc/nS9LKemSN/arHfD+v12tMn
nmsW9Y3GFn661QBvDg0+fLDv1aU25PtTWfd18+DLQ/ADKHKNMXhgrBn7A+nR
yfDvo30dOG8uWasjLo60dcKbR3nD9+68ftsKbz71Gr4MFgV1wFtAI75F/wo2
dMJbKPnsLUb9Xm8Jnf39Nbeg35u/hP78uuo7lN0FP5hq1/Y+Hn52BPtDqLH9
nGs0pBPnhVJOSdmNvMYueGGKF0ZdFfbc5/2H4YVLfewtlfrYi6B5x6wnXoy3
w1tGp39om95sb4Wnx+xFUqh/xqm3XS3wlsv9ZG8FLXptmvi0wYui8aq+aR96
muBFUXRL4pkP6U3wVqL/BnjRpGWuiBxOOQgvRvazFyv7OT+WqgKv2D9Xeuch
FvNUr01uT4gTj8+Lo+KNN+8PHGiGFy/zxF48TY5LUDs8fZ3rSZD+2EsQH1Ov
PH8D8bzXwTNgvmrhrZL5ZW8VnscBeIkSs5eIea2Bl0iGR/HhHzt3w0uSfPaS
UO9+eMnoz+slY16q4SWT3+WzrV/dNfBSMI918FLkfPaMdC4mbejd3Up4RvSz
C55RuX+pch57qZLPXpr8vrOXhuddAS9N+mEvXek3XfpjLwPzuRdehvTHXob4
7GkSs6eh3lJ4JPsxDrLOPkk/vJ/o2gXPC2MPziPleZhknT0T5mknPD1mz4T7
XwbPhHkvh5eJ/kvgZSrvU32dPT1mL1PmiT0z+inA/TEr9ZlRfxE8s1KfHrNn
Ud73Frq4Izn/fcA2eBbMUzE8C/ZvhqfH7Flxfg48q9TLnlXy2bMqnpWm/j2y
KfXZpD/2bJLPnk3xbEp9dvSTC88u/bFnV+qzSz/sedez4a1GfevheWMLPH2d
vdU4zwhPz2fPgXUnPIfiOWQ/ew6phz2HUp9T8fSYPT1mzyn1sueU89hzKffP
pXgumvr/hUvp16XU58a8b4Lnlnz2vPFyeG7x2XNLvexlKf1mKfVlKV6WUp8e
/wdtW/wh
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  ImageSize->{515.0000000000002, Automatic},
  PlotRange->{{0, 60}, {0, 0.2}},
  PlotRangeClipping->True]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10002.},
 CellChangeTimes->{{3.508580842875*^9, 3.50858086753125*^9}, {
  3.508714656546875*^9, 3.508714680796875*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"pdata", "=", 
     RowBox[{"sps", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "1"}], "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"templ", "=", 
     RowBox[{"sps", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tg", "=", 
     RowBox[{"ConstantArray", "[", 
      RowBox[{
       RowBox[{"ConstantArray", "[", 
        RowBox[{"0", ",", 
         RowBox[{"Length", "[", "templ", "]"}]}], "]"}], ",", 
       RowBox[{"Length", "[", "pdata", "]"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"i", "=", "1"}], ";", 
       RowBox[{"k", "=", "1"}], ";", 
       RowBox[{"n", "=", "1"}]}], ",", 
      RowBox[{"i", "\[LessEqual]", 
       RowBox[{"Length", "[", "pdata", "]"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"m", "=", "1"}], ",", 
         RowBox[{"m", "\[LessEqual]", "1"}], ",", 
         RowBox[{"m", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"tg", "[", 
            RowBox[{"[", 
             RowBox[{"n", ",", "k"}], "]"}], "]"}], "=", "1"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"n", "++"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"k", "++"}], ";"}]}], "\[IndentingNewLine]", "]"}], 
    ";"}]}]}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10003.},
 CellChangeTimes->{{3.508581946671875*^9, 3.508581975203125*^9}, {
   3.508582012265625*^9, 3.508582014515625*^9}, 3.508646240125*^9, 
   3.508715233328125*^9, {3.508716461390625*^9, 3.508716463265625*^9}, {
   3.5088841111875*^9, 3.50888412490625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"pg", "=", 
   RowBox[{"classifyDist1", "[", 
    RowBox[{"pdata", ",", "templ", ",", "6.729"}], "]"}]}], ";"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10003.},
 CellChangeTimes->{{3.508581946671875*^9, 3.508581975203125*^9}, {
   3.508582012265625*^9, 3.508582071796875*^9}, {3.508582118421875*^9, 
   3.50858216715625*^9}, {3.508582402953125*^9, 3.50858247465625*^9}, 
   3.508646240125*^9, {3.508715233328125*^9, 3.50871525790625*^9}}],

Cell[BoxData[
 RowBox[{"VF", "[", 
  RowBox[{"tg", ",", "pg"}], "]"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10003.},
 CellChangeTimes->{{3.508582078109375*^9, 3.50858208509375*^9}, 
   3.508715233328125*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.14`", ",", "0.14`", ",", "0.14`"}], "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10003.},
 CellChangeTimes->{
  3.5085820861875*^9, {3.508582120796875*^9, 3.508582168828125*^9}, {
   3.5085824059375*^9, 3.508582476203125*^9}, 3.508715233328125*^9, 
   3.508716497546875*^9, 3.50888431740625*^9, 3.508884361203125*^9}]
}, Closed]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "relation", " ", "of", " ", "number", " ", "of", " ", "features", " ", 
    "and", " ", "EER"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"geteer", "[", 
      RowBox[{"s_", ",", "tg_", ",", "fc_"}], "]"}], ":=", 
     RowBox[{"Module", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"pdata", ",", "templ"}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"pdata", "=", 
         RowBox[{"s", "[", 
          RowBox[{"[", 
           RowBox[{"All", ",", "1"}], "]"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"templ", "=", 
         RowBox[{"s", "[", 
          RowBox[{"[", 
           RowBox[{"All", ",", "2"}], "]"}], "]"}]}], ";", 
        "\[IndentingNewLine]", "\[IndentingNewLine]", 
        RowBox[{"EER", "[", 
         RowBox[{"pdata", ",", "templ", ",", "tg", ",", "fc", ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "10", ",", "0.005"}], "}"}]}], "]"}]}]}], 
      "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tg", "=", 
     RowBox[{"ConstantArray", "[", 
      RowBox[{
       RowBox[{"ConstantArray", "[", 
        RowBox[{"0", ",", 
         RowBox[{"Length", "[", "templ", "]"}]}], "]"}], ",", 
       RowBox[{"Length", "[", "pdata", "]"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"i", "=", "1"}], ";", 
       RowBox[{"k", "=", "1"}], ";", 
       RowBox[{"n", "=", "1"}]}], ",", 
      RowBox[{"i", "\[LessEqual]", 
       RowBox[{"Length", "[", "pdata", "]"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"m", "=", "1"}], ",", 
         RowBox[{"m", "\[LessEqual]", "1"}], ",", 
         RowBox[{"m", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"tg", "[", 
            RowBox[{"[", 
             RowBox[{"n", ",", "k"}], "]"}], "]"}], "=", "1"}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"n", "++"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"k", "++"}], ";"}]}], "\[IndentingNewLine]", "]"}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5088250114375*^9, 3.508825074203125*^9}, {
   3.50883141278125*^9, 3.508831454796875*^9}, 3.508831485625*^9, {
   3.508834868734375*^9, 3.50883492403125*^9}, {3.508838791140625*^9, 
   3.50883879328125*^9}, {3.508838833890625*^9, 3.50883903178125*^9}, 
   3.508839169421875*^9, {3.508839201875*^9, 3.508839221609375*^9}, {
   3.508839397890625*^9, 3.5088394080625*^9}, {3.508839526*^9, 
   3.50883953828125*^9}, {3.508842065515625*^9, 3.508842111515625*^9}, {
   3.50888402603125*^9, 3.5088840274375*^9}, 3.508884140015625*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"n", "=", "1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ps", "=", 
   RowBox[{"Partition", "[", 
    RowBox[{
     RowBox[{"tpc", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", 
        RowBox[{"1", ";;", "n"}]}], "]"}], "]"}], ",", "2"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Dimensions", "[", "ps", "]"}]}], "Input",
 CellChangeTimes->{{3.5088250114375*^9, 3.508825074203125*^9}, {
   3.50883141278125*^9, 3.508831454796875*^9}, 3.508831485625*^9, {
   3.508834868734375*^9, 3.50883492403125*^9}, {3.508838791140625*^9, 
   3.50883879328125*^9}, {3.508838833890625*^9, 3.508839037296875*^9}, {
   3.508839438140625*^9, 3.508839442078125*^9}, {3.508839575328125*^9, 
   3.508839576046875*^9}, 3.508839785109375*^9, 3.50884203296875*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"50", ",", "2", ",", "1"}], "}"}]], "Output",
 CellChangeTimes->{
  3.508839442796875*^9, 3.508839577453125*^9, 3.50884204034375*^9, {
   3.508842102484375*^9, 3.50884211675*^9}, 3.508884374890625*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"e", "=", 
    RowBox[{
     RowBox[{"Reap", "[", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"i", "=", "1"}], ",", 
        RowBox[{"i", "\[LessEqual]", 
         RowBox[{"Length", "[", 
          RowBox[{"tpc", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "]"}]}], ",", 
        RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"Sow", "[", 
          RowBox[{"Mean", "[", 
           RowBox[{
            RowBox[{"geteer", "[", 
             RowBox[{
              RowBox[{"Partition", "[", 
               RowBox[{
                RowBox[{"tpc", "[", 
                 RowBox[{"[", 
                  RowBox[{";;", ",", 
                   RowBox[{";;", "i"}]}], "]"}], "]"}], ",", "2"}], "]"}], 
              ",", "tg", ",", "classifyDist1"}], "]"}], "[", 
            RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"Print", "[", "i", "]"}], ";"}]}], "]"}], "]"}], "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ";"}], "*)"}]], "Input",
 CellChangeTimes->{{3.508839044078125*^9, 3.508839060484375*^9}, 
   3.508839190734375*^9, {3.508839234421875*^9, 3.508839241296875*^9}, {
   3.508839530109375*^9, 3.508839531953125*^9}, {3.508839702578125*^9, 
   3.50883977965625*^9}, 3.50883981765625*^9, {3.50884190090625*^9, 
   3.50884195915625*^9}, {3.508842220625*^9, 3.50884224371875*^9}, {
   3.5088423219375*^9, 3.508842322515625*^9}, {3.508842905609375*^9, 
   3.5088429259375*^9}, 3.508883650890625*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"ee", "=", 
   RowBox[{"Import", "[", "\"\<Cache_edat.mx\>\"", "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.508846834625*^9, 3.50884688190625*^9}, {
  3.5088795541875*^9, 3.508879568921875*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Dimensions", "[", "ee", "]"}], "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.5088795749375*^9, 3.5088795778125*^9}, {
   3.508879611828125*^9, 3.508879612328125*^9}, 3.508884390515625*^9}],

Cell[BoxData[
 RowBox[{"{", "99", "}"}]], "Output",
 CellChangeTimes->{3.508879578609375*^9, 3.508879613625*^9, 
  3.50888439221875*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"100", "*", "ee"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"5", ",", "20"}], "}"}]}], ",", 
   RowBox[{"Mesh", "\[Rule]", "Full"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5088795848125*^9, 3.508879588546875*^9}, {
  3.508879739484375*^9, 3.5088798006875*^9}, {3.50887984334375*^9, 
  3.5088798475625*^9}, {3.508902755859375*^9, 3.508902764578125*^9}}],

Cell[BoxData[
 GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJzt1blKA1EYxfGjIIiNFilEBDUEFQnucddcjftu4pJaCdgJ9sIFBStrH8HG
h7hVaitrG/s0NlZmuReZ/zPcgeHjN2fODMMU39D1XbHSLqlaPxuzddTyuc/n
6u9H1vgLfmb87DSPr6XafVdwD/KUebh97375GfHu9XPUz3540LQ1Z3hfGs6Y
p+b7goeRh+dM+jkGZ+FxmHMKnoZn4FlTuGl8cXAOnoPn4QV40XRMvH19p4KX
4GV4xdR/VnmgL3gV+Rqcx/0maQtrPdm3sDaQwyogh7WJHNYWcljbyGHtIIe1
ixzWHnJY+8hhHSCHdYgc1hFyWMfIYZ0gh3WatIEt7GCdoQ9b2MEqog9b2MEq
oQ9b2ME6Rx+2sIN1gT5sYQfrEn3Ywg7WFfqwhR2sMvqwhR0c9wUc90XcF3Ff
xH0R90Wr/+8/P5m0Hg==
   "], {{{}, {}, 
     {Hue[0.67, 0.6, 0.6], 
      LineBox[{100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 
       112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 
       126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 
       140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 
       154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 
       168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 
       182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 
       196, 197, 198}]}}, {
     {Hue[0.67, 0.6, 0.6], 
      PointBox[{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18,
        19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 
       36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,
        54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 
       71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88,
        89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99}]}, {}, {}}}],
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 6.},
  ImageSize->{514.1666666666666, Automatic},
  PlotRange->{{0, 99.}, {5, 20}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], Automatic}]], "Output",
 CellChangeTimes->{
  3.5088795895625*^9, {3.50887974125*^9, 3.5088797645*^9}, 
   3.50887980146875*^9, 3.508879847796875*^9, 3.508884394609375*^9, 
   3.5089027655*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
  "select", " ", "more", " ", "effective", " ", "pcs", " ", "use", " ", 
   "MDA"}], "*)"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10004.},
 CellChangeTimes->{{3.508646248828125*^9, 3.50864626675*^9}, 
   3.508715281640625*^9, {3.508884184328125*^9, 3.508884187796875*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"spc", "=", 
   RowBox[{"Partition", "[", 
    RowBox[{"tpc", ",", "2"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Dimensions", "[", "%", "]"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10004.},
 CellChangeTimes->{
  3.508646282828125*^9, {3.508646314515625*^9, 3.5086463309375*^9}, 
   3.508715281640625*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"50", ",", "2", ",", "99"}], "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10004.},
 CellChangeTimes->{3.508649823546875*^9, 3.508715281640625*^9, 
  3.50871531165625*^9, 3.508884400046875*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"scatter", " ", "matrix", " ", "of", " ", "all", " ", "pcs"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Sw", "=", 
     RowBox[{"Tr", "[", 
      RowBox[{"Total", "[", 
       RowBox[{"Map", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Covariance", "[", "#", "]"}], "&"}], ",", "spc"}], "]"}], 
       "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Dimensions", "[", "%", "]"}], "\[IndentingNewLine]", "m"}], 
     "=", 
     RowBox[{"Mean", "[", "tpc", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"mm", "=", 
     RowBox[{"Map", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Mean", "[", "#", "]"}], "&"}], ",", "spc"}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Sb", "=", 
     RowBox[{"Total", "[", 
      RowBox[{"Map", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"#", "-", "m"}], ")"}], ".", 
          RowBox[{"(", 
           RowBox[{"#", "-", "m"}], ")"}]}], "&"}], ",", "mm"}], "]"}], 
      "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Dimensions", "[", "%", "]"}]}], "*)"}]}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.508646342171875*^9, 3.50864637528125*^9}, {
   3.508646412140625*^9, 3.5086464904375*^9}, {3.508646535453125*^9, 
   3.508646583203125*^9}, {3.508646660984375*^9, 3.508646665875*^9}, {
   3.5086474775*^9, 3.508647479890625*^9}, {3.508647563984375*^9, 
   3.508647580109375*^9}, {3.508649831921875*^9, 3.508649847796875*^9}, 
   3.5086513170625*^9, 3.5086513493125*^9}],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.508649833203125*^9, 3.508649849015625*^9}, 
   3.5086513493125*^9}],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.508649833203125*^9, 3.50864984903125*^9}, 
   3.5086513493125*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"spc2", "=", 
    RowBox[{"Partition", "[", 
     RowBox[{
      RowBox[{"tpc", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "1"}], "]"}], "]"}], ",", "2"}], "]"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"Dimensions", "[", "%", "]"}]}], "*)"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.508647589*^9, 3.508647603265625*^9}, {
  3.508647648171875*^9, 3.508647648578125*^9}, {3.508651329703125*^9, 
  3.508651337078125*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"50", ",", "2"}], "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{3.508647651109375*^9, 3.508651337078125*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"Sw", "=", 
    RowBox[{"Total", "[", 
     RowBox[{"Map", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Covariance", "[", "#", "]"}], "&"}], ",", "spc2"}], "]"}], 
     "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"m", "=", 
    RowBox[{"Mean", "[", 
     RowBox[{"tpc", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "1"}], "]"}], "]"}], "]"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"Flatten", "[", 
    RowBox[{"spc2", ",", "1"}], "]"}], ";", "\[IndentingNewLine]", 
   RowBox[{"mm", "=", 
    RowBox[{"Map", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"Mean", "[", "#", "]"}], "&"}], ",", "spc2"}], "]"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"Mean", "[", "mm", "]"}], ";", "\[IndentingNewLine]", 
   RowBox[{"Sb", "=", 
    RowBox[{"Total", "[", 
     RowBox[{"Map", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"#", "-", "m"}], ")"}], "*", 
         RowBox[{"(", 
          RowBox[{"#", "-", "m"}], ")"}]}], "&"}], ",", "mm"}], "]"}], 
     "]"}]}], ";"}], "*)"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.508647673640625*^9, 3.50864773103125*^9}, {
   3.508648724484375*^9, 3.50864872471875*^9}, {3.5086488176875*^9, 
   3.508648817859375*^9}, {3.508649000515625*^9, 3.508649066125*^9}, {
   3.508649381703125*^9, 3.508649387546875*^9}, 3.50864942815625*^9, 
   3.508651337078125*^9}]
}, Closed]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "this", " ", "will", " ", "not", " ", "change", " ", "even", " ", 
    "standardized", " ", "PC"}], "*)"}], " ", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"sj1", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"J1", "[", 
       RowBox[{"Partition", "[", 
        RowBox[{
         RowBox[{"tpc", "[", 
          RowBox[{"[", 
           RowBox[{"All", ",", "i"}], "]"}], "]"}], ",", "2"}], "]"}], "]"}], 
      ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"Length", "[", 
         RowBox[{"tpc", "[", 
          RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "}"}]}], "]"}]}], 
   ";"}]}]], "Input",
 CellChangeTimes->{{3.508651419421875*^9, 3.508651476765625*^9}, {
  3.51005126821875*^9, 3.51005129715625*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"sj1", ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"Mesh", "\[Rule]", "Full"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"0", ",", "25"}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.508649395796875*^9, 3.508649411296875*^9}, {
  3.5086494995*^9, 3.508649609046875*^9}, {3.50865143465625*^9, 
  3.508651436453125*^9}}],

Cell[BoxData[
 GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJztlftT1GUUxlekAGkCFYSUEnJBIpCLgIDgPqzCwiLLsuzCEkghZdhEIBdF
zWkD5CaVgVTTGGCSg44oBkYkjAsIFPeLcStwuQUIGCUYEmPZcN4f+w92Z76z
8/k+53nf97zn7B6zg7GStzQ4HE7j0+e/79XPPC87PSty7nVn0AsEHUN0Ur8Z
sTZuKrn200YbifUxmWEe7VW8htgAfV9bx+SnzfNW2Riq4vkP85ruE5vA4tLc
+gnrMWJTbK2O6wodHiR+GZbfi+tuVGrSelzoXiq/rnFxgXQLdEeffLaHN05s
iRWPX4X105PEVsj8UTywPW2F2BpTRWMb4wKmiXdA5XTCMnuJxdshdlv1CYfR
KWJ79Od4Rql2M90Bh6URb7+xY4Z4JxLcLZKlx38ndoT5+IJ+VPoosRNiHxul
TX6xSOyMNXFaX3nHsPV2Qa9RFjlQx/JxwR4t3aPJZxm7ouGfSNnmFBbvhsbC
+Z/K+Gz93RjLFL5SFPGA2B3yjGMPs1pniT1Q3XurmafH9t+D3FhDt4PPdRPz
EOouqqvuHCAGntnWUVkgGlplBTBxpmdOR95HuieMz94tsfWg/RWe6DWc/6hp
5DfS+Sgz4JrYGt0jnQ/deO3LL33LzrsXOtrmPgmazL8XVg4GXNl6iufsg2NL
vNbza1Wk74OPdlWBUeoI6V64X29jN2NK51d4oa3OZUQ000O6Nw4kCVVv2lC8
whuXxyXZxT/0ki7A5rGyi5/YKUkXYAP3r1r/SdZvPvCUCtaGDrP8fXCu7uq5
CsO7pPvC9ANB3OFupvsiYvqKUJrLdCH4hvKctIZfSBeC57qYHnGT3bcfOpSz
o4O2LaT7wTFX1HYyj+n74ffxsnPiNaqHYj+S7uVGtLnSfXD8UaN9xi7vVdpP
4Y9w+VjpbTfmF+FBfGVa4mnmFyFu1HvnqVJ2PwE4bvBemJkm8weAcyHq760F
E6SL8bOTUbNVLtUbYkSZFPtWlLB4Mdysa1rydZpXWSnGSOqphc4XWX0Dkf/Y
0rxnE/0eEYilwpXXBOOt5A+E1rK9ys6W8lE+5Vu90w09rP8k4Kdygq99R/0A
CW4Uljl85sn6QQLd/hcsEmvpfEoJHrY2WZjwmT8Io+KMA9w8qj+C0F2eobPO
gOIVQTiUqJq6U0LxyiA8CWuqkZWy/pLC/0jL+T4vlr8UKdn8xS2PusgvxTuf
fh7u4kT/V0ophrNEzn/Ws/rLELUs1z+qy/wyvM+9vSlikPWjDBM+V85/k0z5
KWWwzBMXPslh/mAMRmpkh01RvRCM7a4VHlW8dvIHY1d2iH76eAf5g+HAGWpd
t6WW/CFwNO5KeHeM8kMIirsN0mO+pPUUIeirmmz341M9lCF4VCSxnr1aR345
/riTueFCJ9UXcugNHKoeCmgkvxw24cPXjyTRfSjlaE9Zmha4lZNfPS/U80I9
L9TzQj0vVv3qefF/8+Jfo3mumw==
   "], {{{}, {}, 
     {Hue[0.67, 0.6, 0.6], 
      LineBox[{100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 
       112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 
       126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 
       140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 
       154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 
       168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 
       182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 
       196, 197, 198}]}}, {
     {Hue[0.67, 0.6, 0.6], 
      PointBox[{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18,
        19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 
       36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,
        54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 
       71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88,
        89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99}]}, {}, {}}}],
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  ImageSize->{451.6666666666664, Automatic},
  PlotRange->{{0, 99.}, {0, 25}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], Automatic}]], "Output",
 CellChangeTimes->{
  3.508649412125*^9, {3.5086495526875*^9, 3.508649609984375*^9}, 
   3.508649877796875*^9, 3.50865000915625*^9, 3.508650549796875*^9, 
   3.50865144025*^9, 3.508715336734375*^9, 3.50888441178125*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"features", " ", "which", " ", "have", " ", "j1"}], ">", "1"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"ind", "=", 
    RowBox[{"Reverse", "[", 
     RowBox[{"Ordering", "[", "sj1", "]"}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Length", "[", "ind", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.5086514875*^9, 3.508651550375*^9}, {
   3.508651830890625*^9, 3.5086518785625*^9}, {3.508651970296875*^9, 
   3.5086519936875*^9}, 3.5087163040625*^9, {3.508884425625*^9, 
   3.508884433921875*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1", ",", "2", ",", "3", ",", "9", ",", "4", ",", "13", ",", "21", ",", 
   "29", ",", "23", ",", "10", ",", "24", ",", "5", ",", "19", ",", "27", ",",
    "28", ",", "18", ",", "6", ",", "14", ",", "16", ",", "15", ",", "12", 
   ",", "25", ",", "22", ",", "17", ",", "35", ",", "63", ",", "11", ",", 
   "68", ",", "82", ",", "7", ",", "37", ",", "20", ",", "26", ",", "34", ",",
    "67", ",", "38", ",", "85", ",", "79", ",", "76", ",", "41", ",", "44", 
   ",", "73", ",", "57", ",", "70", ",", "40", ",", "72", ",", "86", ",", 
   "36", ",", "39", ",", "48", ",", "50", ",", "32", ",", "52", ",", "47", 
   ",", "8", ",", "92", ",", "71", ",", "31", ",", "56", ",", "78", ",", "60",
    ",", "75", ",", "74", ",", "64", ",", "33", ",", "80", ",", "77", ",", 
   "84", ",", "45", ",", "83", ",", "65", ",", "51", ",", "87", ",", "58", 
   ",", "49", ",", "62", ",", "43", ",", "88", ",", "61", ",", "93", ",", 
   "53", ",", "55", ",", "59", ",", "42", ",", "30", ",", "81", ",", "98", 
   ",", "90", ",", "89", ",", "69", ",", "94", ",", "54", ",", "96", ",", 
   "66", ",", "97", ",", "95", ",", "91", ",", "46", ",", "99"}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.508651985625*^9, 3.5086519946875*^9}, 
   3.508716305546875*^9, 3.508884436109375*^9}],

Cell[BoxData["99"], "Output",
 CellChangeTimes->{{3.508651985625*^9, 3.5086519946875*^9}, 
   3.508716305546875*^9, 3.508884436109375*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"stpc", "=", 
    RowBox[{"tpc", "[", 
     RowBox[{"[", 
      RowBox[{"All", ",", "ind"}], "]"}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.50865189671875*^9, 3.508651920859375*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"J1", "[", 
  RowBox[{"Partition", "[", 
   RowBox[{"stpc", ",", "2"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.508649770328125*^9, 3.508649788203125*^9}, {
  3.508650015125*^9, 3.508650015390625*^9}, {3.508651257515625*^9, 
  3.50865126153125*^9}, {3.508651951859375*^9, 3.50865195246875*^9}, {
  3.5088844826875*^9, 3.5088845299375*^9}}],

Cell[BoxData["9.749950252866475`"], "Output",
 CellChangeTimes->{{3.508884510640625*^9, 3.508884530671875*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
  "what", " ", "if", " ", "pcs", " ", "are", " ", "ordered", " ", "by", " ", 
   "MDA"}], "*)"}]], "Input",
 CellChangeTimes->{{3.508884926*^9, 3.508884944890625*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"t1", "=", 
    RowBox[{"AbsoluteTime", "[", "]"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"e", "=", 
    RowBox[{
     RowBox[{"Reap", "[", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"i", "=", "1"}], ",", 
        RowBox[{"i", "\[LessEqual]", 
         RowBox[{"Length", "[", 
          RowBox[{"stpc", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "]"}]}], ",", 
        RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"Sow", "[", 
          RowBox[{"Mean", "[", 
           RowBox[{
            RowBox[{"geteer", "[", 
             RowBox[{
              RowBox[{"Partition", "[", 
               RowBox[{
                RowBox[{"stpc", "[", 
                 RowBox[{"[", 
                  RowBox[{";;", ",", 
                   RowBox[{";;", "i"}]}], "]"}], "]"}], ",", "2"}], "]"}], 
              ",", "tg", ",", "classifyDist1"}], "]"}], "[", 
            RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"Print", "[", "i", "]"}], ";"}]}], "]"}], "]"}], "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
   RowBox[{"t2", "=", 
    RowBox[{
     RowBox[{"AbsoluteTime", "[", "]"}], "-", "t1"}]}]}], "*)"}]], "Input",
 CellChangeTimes->{{3.508884898890625*^9, 3.5088849054375*^9}, {
   3.50888497290625*^9, 3.508884987390625*^9}, {3.508885044234375*^9, 
   3.5088850594375*^9}, 3.508905109265625*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"e", "=", 
   RowBox[{"Import", "[", "\"\<Cache_ePcamda.mx\>\"", "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.508905113*^9, 3.50890512659375*^9}}],

Cell[CellGroupData[{

Cell[BoxData["1"], "Print",
 CellChangeTimes->{3.508885083125*^9}],

Cell[BoxData["2"], "Print",
 CellChangeTimes->{3.5088851036875*^9}],

Cell[BoxData["3"], "Print",
 CellChangeTimes->{3.50888512646875*^9}],

Cell[BoxData["4"], "Print",
 CellChangeTimes->{3.508885149765625*^9}],

Cell[BoxData["5"], "Print",
 CellChangeTimes->{3.508885174203125*^9}],

Cell[BoxData["6"], "Print",
 CellChangeTimes->{3.50888519903125*^9}],

Cell[BoxData["7"], "Print",
 CellChangeTimes->{3.508885224359375*^9}],

Cell[BoxData["8"], "Print",
 CellChangeTimes->{3.50888524915625*^9}],

Cell[BoxData["9"], "Print",
 CellChangeTimes->{3.5088852744375*^9}],

Cell[BoxData["10"], "Print",
 CellChangeTimes->{3.508885300109375*^9}],

Cell[BoxData["11"], "Print",
 CellChangeTimes->{3.508885325828125*^9}],

Cell[BoxData["12"], "Print",
 CellChangeTimes->{3.508885354484375*^9}],

Cell[BoxData["13"], "Print",
 CellChangeTimes->{3.508885383375*^9}],

Cell[BoxData["14"], "Print",
 CellChangeTimes->{3.50888541209375*^9}],

Cell[BoxData["15"], "Print",
 CellChangeTimes->{3.508885440703125*^9}],

Cell[BoxData["16"], "Print",
 CellChangeTimes->{3.50888547265625*^9}],

Cell[BoxData["17"], "Print",
 CellChangeTimes->{3.50888550546875*^9}],

Cell[BoxData["18"], "Print",
 CellChangeTimes->{3.508885535484375*^9}],

Cell[BoxData["19"], "Print",
 CellChangeTimes->{3.508885565515625*^9}],

Cell[BoxData["20"], "Print",
 CellChangeTimes->{3.50888559584375*^9}],

Cell[BoxData["21"], "Print",
 CellChangeTimes->{3.508885626375*^9}],

Cell[BoxData["22"], "Print",
 CellChangeTimes->{3.508885657109375*^9}],

Cell[BoxData["23"], "Print",
 CellChangeTimes->{3.508885688046875*^9}],

Cell[BoxData["24"], "Print",
 CellChangeTimes->{3.508885718609375*^9}],

Cell[BoxData["25"], "Print",
 CellChangeTimes->{3.508885749421875*^9}],

Cell[BoxData["26"], "Print",
 CellChangeTimes->{3.50888578040625*^9}],

Cell[BoxData["27"], "Print",
 CellChangeTimes->{3.50888581115625*^9}],

Cell[BoxData["28"], "Print",
 CellChangeTimes->{3.508885842296875*^9}],

Cell[BoxData["29"], "Print",
 CellChangeTimes->{3.508885873234375*^9}],

Cell[BoxData["30"], "Print",
 CellChangeTimes->{3.508885904203125*^9}],

Cell[BoxData["31"], "Print",
 CellChangeTimes->{3.508885935703125*^9}],

Cell[BoxData["32"], "Print",
 CellChangeTimes->{3.50888596696875*^9}],

Cell[BoxData["33"], "Print",
 CellChangeTimes->{3.50888599815625*^9}],

Cell[BoxData["34"], "Print",
 CellChangeTimes->{3.508886029203125*^9}],

Cell[BoxData["35"], "Print",
 CellChangeTimes->{3.508886060390625*^9}],

Cell[BoxData["36"], "Print",
 CellChangeTimes->{3.508886091375*^9}],

Cell[BoxData["37"], "Print",
 CellChangeTimes->{3.508886122296875*^9}],

Cell[BoxData["38"], "Print",
 CellChangeTimes->{3.5088861531875*^9}],

Cell[BoxData["39"], "Print",
 CellChangeTimes->{3.50888618409375*^9}],

Cell[BoxData["40"], "Print",
 CellChangeTimes->{3.508886217953125*^9}],

Cell[BoxData["41"], "Print",
 CellChangeTimes->{3.50888625184375*^9}],

Cell[BoxData["42"], "Print",
 CellChangeTimes->{3.508886285671875*^9}],

Cell[BoxData["43"], "Print",
 CellChangeTimes->{3.5088863195*^9}],

Cell[BoxData["44"], "Print",
 CellChangeTimes->{3.5088863505625*^9}],

Cell[BoxData["45"], "Print",
 CellChangeTimes->{3.508886381671875*^9}],

Cell[BoxData["46"], "Print",
 CellChangeTimes->{3.508886412671875*^9}],

Cell[BoxData["47"], "Print",
 CellChangeTimes->{3.508886443703125*^9}],

Cell[BoxData["48"], "Print",
 CellChangeTimes->{3.508886478640625*^9}],

Cell[BoxData["49"], "Print",
 CellChangeTimes->{3.508886512875*^9}],

Cell[BoxData["50"], "Print",
 CellChangeTimes->{3.50888654771875*^9}],

Cell[BoxData["51"], "Print",
 CellChangeTimes->{3.508886583328125*^9}],

Cell[BoxData["52"], "Print",
 CellChangeTimes->{3.508886615859375*^9}],

Cell[BoxData["53"], "Print",
 CellChangeTimes->{3.508886648046875*^9}],

Cell[BoxData["54"], "Print",
 CellChangeTimes->{3.508886680515625*^9}],

Cell[BoxData["55"], "Print",
 CellChangeTimes->{3.50888671759375*^9}],

Cell[BoxData["56"], "Print",
 CellChangeTimes->{3.508886755*^9}],

Cell[BoxData["57"], "Print",
 CellChangeTimes->{3.5088867928125*^9}],

Cell[BoxData["58"], "Print",
 CellChangeTimes->{3.508886830078125*^9}],

Cell[BoxData["59"], "Print",
 CellChangeTimes->{3.50888686734375*^9}],

Cell[BoxData["60"], "Print",
 CellChangeTimes->{3.508886904828125*^9}],

Cell[BoxData["61"], "Print",
 CellChangeTimes->{3.5088869421875*^9}],

Cell[BoxData["62"], "Print",
 CellChangeTimes->{3.50888698059375*^9}],

Cell[BoxData["63"], "Print",
 CellChangeTimes->{3.508887017796875*^9}],

Cell[BoxData["64"], "Print",
 CellChangeTimes->{3.508887058875*^9}],

Cell[BoxData["65"], "Print",
 CellChangeTimes->{3.508887100125*^9}],

Cell[BoxData["66"], "Print",
 CellChangeTimes->{3.5088871411875*^9}],

Cell[BoxData["67"], "Print",
 CellChangeTimes->{3.508887182515625*^9}],

Cell[BoxData["68"], "Print",
 CellChangeTimes->{3.508887224046875*^9}],

Cell[BoxData["69"], "Print",
 CellChangeTimes->{3.508887265796875*^9}],

Cell[BoxData["70"], "Print",
 CellChangeTimes->{3.508887307*^9}],

Cell[BoxData["71"], "Print",
 CellChangeTimes->{3.508887348609375*^9}],

Cell[BoxData["72"], "Print",
 CellChangeTimes->{3.508887386484375*^9}],

Cell[BoxData["73"], "Print",
 CellChangeTimes->{3.50888742428125*^9}],

Cell[BoxData["74"], "Print",
 CellChangeTimes->{3.508887462*^9}],

Cell[BoxData["75"], "Print",
 CellChangeTimes->{3.508887499796875*^9}],

Cell[BoxData["76"], "Print",
 CellChangeTimes->{3.5088875375625*^9}],

Cell[BoxData["77"], "Print",
 CellChangeTimes->{3.5088875753125*^9}],

Cell[BoxData["78"], "Print",
 CellChangeTimes->{3.508887613046875*^9}],

Cell[BoxData["79"], "Print",
 CellChangeTimes->{3.50888765090625*^9}],

Cell[BoxData["80"], "Print",
 CellChangeTimes->{3.508887688890625*^9}],

Cell[BoxData["81"], "Print",
 CellChangeTimes->{3.508887726640625*^9}],

Cell[BoxData["82"], "Print",
 CellChangeTimes->{3.50888776446875*^9}],

Cell[BoxData["83"], "Print",
 CellChangeTimes->{3.508887802265625*^9}],

Cell[BoxData["84"], "Print",
 CellChangeTimes->{3.508887840078125*^9}],

Cell[BoxData["85"], "Print",
 CellChangeTimes->{3.5088878779375*^9}],

Cell[BoxData["86"], "Print",
 CellChangeTimes->{3.508887915890625*^9}],

Cell[BoxData["87"], "Print",
 CellChangeTimes->{3.50888795403125*^9}],

Cell[BoxData["88"], "Print",
 CellChangeTimes->{3.508887992375*^9}],

Cell[BoxData["89"], "Print",
 CellChangeTimes->{3.508888030140625*^9}],

Cell[BoxData["90"], "Print",
 CellChangeTimes->{3.50888806828125*^9}],

Cell[BoxData["91"], "Print",
 CellChangeTimes->{3.50888810646875*^9}],

Cell[BoxData["92"], "Print",
 CellChangeTimes->{3.508888144546875*^9}],

Cell[BoxData["93"], "Print",
 CellChangeTimes->{3.508888182765625*^9}],

Cell[BoxData["94"], "Print",
 CellChangeTimes->{3.508888219796875*^9}],

Cell[BoxData["95"], "Print",
 CellChangeTimes->{3.5088882566875*^9}],

Cell[BoxData["96"], "Print",
 CellChangeTimes->{3.508888293546875*^9}],

Cell[BoxData["97"], "Print",
 CellChangeTimes->{3.508888330328125*^9}],

Cell[BoxData["98"], "Print",
 CellChangeTimes->{3.508888367171875*^9}],

Cell[BoxData["99"], "Print",
 CellChangeTimes->{3.508888404828125*^9}]
}, Closed]],

Cell[BoxData["3337.296875`10.974939835283967"], "Output",
 CellChangeTimes->{3.50888840484375*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"Export", "[", 
    RowBox[{"\"\<Cache_ePcamda.mx\>\"", ",", "e"}], "]"}], ";"}], 
  "*)"}]], "Input",
 CellChangeTimes->{{3.508903780609375*^9, 3.508903798109375*^9}, 
   3.508905105875*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"The", " ", "result", " ", "seems", " ", "much", " ", "better"}], 
   "*)"}], "\[IndentingNewLine]"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.508902841171875*^9, 3.508902848890625*^9}, {
   3.508903126890625*^9, 3.508903142109375*^9}, 3.508905141515625*^9, 
   3.509598161375*^9}],

Cell[BoxData[{
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"100", "*", "e"}], ",", 
   RowBox[{"Joined", "\[Rule]", "True"}], ",", 
   RowBox[{"Mesh", "\[Rule]", "Full"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"9.5", ",", "16.5"}], "}"}]}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"e", "[", 
  RowBox[{"[", 
   RowBox[{";;", "11"}], "]"}], "]"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.5088892924375*^9, 3.508889300015625*^9}, {
   3.508902711265625*^9, 3.5089027204375*^9}, {3.508902809984375*^9, 
   3.508902839515625*^9}, {3.508903121046875*^9, 3.50890312340625*^9}, {
   3.508905143953125*^9, 3.508905222640625*^9}, 3.509598161375*^9}],

Cell[BoxData[
 GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJzt1blKQ0EUh/E/ghBstEiRwkJFVCS4x10zcde4Jm61INgJ9sIBCytrH8HG
hzhVaqvUNvZpbKw0eFPM9wxz4TL88s25VeAM3tzXbrskNf7e9vn/tMql5nPj
56MYsh+yczg7c+Hptd566Om4D87jfgHux/0BeAjmOQKP4RyHi/AEPAVPh8e7
996X79HMM+iz8BxcgufxvQV4EfeX4GV4BV6F1+B1uAyH2Aargg5rAx3WJjqs
LXRY2+iwdtBh7aLD2kOHtY8O6wAdVhW9Gron3z6/8p3/3WFsg3WEDusYHdYJ
OqzT2AE22GGdYR422GHVMA8b7LDqmIcNdljnmIcNdlgXmIcNdliXmIcNdlhX
mIcNdljXmIcNdjjti7Qvop72RdzTvkj7Iu2Ltn8Bb9USVA==
   "], {{{}, {}, 
     {Hue[0.67, 0.6, 0.6], 
      LineBox[{100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 
       112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 
       126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 
       140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 
       154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 
       168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 
       182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 
       196, 197, 198}]}}, {
     {Hue[0.67, 0.6, 0.6], 
      PointBox[{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18,
        19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 
       36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,
        54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 
       71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88,
        89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99}]}, {}, {}}}],
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 10.},
  ImageSize->{602.5000000000002, Automatic},
  PlotRange->{{0, 99.}, {9.5, 16.5}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], Automatic}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{
  3.508905145640625*^9, {3.508905177234375*^9, 3.508905224140625*^9}, 
   3.509598161375*^9}]
}, Closed]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"use", " ", "50", " ", 
     RowBox[{"persons", "'"}], " ", "one", " ", "measure", " ", "as", " ", 
     "a", " ", "sample"}], ",", " ", 
    RowBox[{"how", " ", "to", " ", "compare", " ", "same", " ", 
     RowBox[{"person", "'"}], "s", " ", 
     RowBox[{"measure", "?"}]}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "leave", " ", "it", " ", "to", " ", "further", " ", "research", " ", 
    "when", " ", "there", " ", "are", " ", "enough", " ", "data"}], 
   "*)"}]}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 CellChangeTimes->{{3.508902821203125*^9, 3.508902831*^9}, 
   3.50890522959375*^9, {3.508905328546875*^9, 3.50890533646875*^9}, {
   3.509167617078125*^9, 3.509167637140625*^9}, {3.509597633984375*^9, 
   3.50959765875*^9}, {3.509597851890625*^9, 3.50959785278125*^9}, {
   3.509597923078125*^9, 3.509597923515625*^9}, {3.5095988395*^9, 
   3.5095988650625*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"fft1", "=", 
   RowBox[{"Flatten", "[", 
    RowBox[{
     RowBox[{"data", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "1"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", "1", "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"2", ",", "3"}], "}"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fft2", "=", 
   RowBox[{"Flatten", "[", 
    RowBox[{
     RowBox[{"data", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "2"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", "1", "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"2", ",", "3"}], "}"}]}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{"Dimensions", "[", "fft1", "]"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.50916765053125*^9, 3.509167660421875*^9}, {
   3.509167706390625*^9, 3.50916772003125*^9}, 3.509598835375*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"50", ",", "276"}], "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{3.5091676635*^9, 3.509167722953125*^9, 
  3.509168614703125*^9, 3.509596499*^9, 3.509598835375*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"sample", "=", "fft1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"test", "=", "fft2"}], ";"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.509596522921875*^9, 3.50959652728125*^9}, {
   3.509597173546875*^9, 3.5095971781875*^9}, 3.509598835375*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"pc1", "=", 
   RowBox[{"PrincipalComponents", " ", "[", "sample", "]"}]}], ";"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{3.509596557453125*^9, 3.509598835375*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"ssam", "=", 
   RowBox[{"StdP", "[", "sample", "]"}]}], ";"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.509596882515625*^9, 3.509596892359375*^9}, 
   3.509598835375*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"tt", ",", "v", ",", "mpc1"}], "}"}], "=", 
   RowBox[{"mPCA", "[", "sample", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Dimensions", "[", "%", "]"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.50916780296875*^9, 3.509167833828125*^9}, {
   3.509168645890625*^9, 3.509168646984375*^9}, {3.50959674975*^9, 
   3.509596760109375*^9}, 3.509598835375*^9}],

Cell[BoxData[
 RowBox[{"Dimensions", "[", "tt", "]"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.509597890015625*^9, 3.509597890203125*^9}, 
   3.509598835375*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"276", ",", "276"}], "}"}]], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{3.509597890890625*^9, 3.509598835375*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Chop", "[", 
   RowBox[{"ssam", ".", 
    RowBox[{"tt", "\[Transpose]"}]}], "]"}], "\[Equal]", "mpc1"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.509596898546875*^9, 3.50959694090625*^9}, 
   3.509598835375*^9}],

Cell[BoxData["True"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{
  3.509168234609375*^9, 3.509168619734375*^9, 3.50916864978125*^9, 
   3.509596531265625*^9, 3.50959676109375*^9, {3.50959690934375*^9, 
   3.5095969415625*^9}, 3.509598835375*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Chop", "[", 
    RowBox[{"mpc1", "-", 
     RowBox[{"ssam", ".", 
      RowBox[{"tt", "\[Transpose]"}]}]}], "]"}], "[", 
   RowBox[{"[", 
    RowBox[{"All", ",", 
     RowBox[{";;", "100"}]}], "]"}], "]"}], ";"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.509167849375*^9, 3.50916787890625*^9}, 
   3.5095965686875*^9, {3.50959692634375*^9, 3.509596929375*^9}, 
   3.509597108984375*^9, 3.509598835375*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"pmm", "=", 
    RowBox[{"Mean", "[", "sample", "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"stest", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"test", "[", 
       RowBox[{"[", "i", "]"}], "]"}], "-", "pmm"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"Length", "[", "test", "]"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ptest", "=", 
   RowBox[{"stest", ".", 
    RowBox[{"tt", "\[Transpose]"}]}]}], ";"}]}], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.50959714509375*^9, 3.5095972118125*^9}, {
   3.509597274546875*^9, 3.50959738021875*^9}, 3.509598835375*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"indd", "=", 
   RowBox[{"{", 
    RowBox[{
    "1", ",", "2", ",", "3", ",", "9", ",", "4", ",", "13", ",", "21", ",", 
     "29", ",", "23", ",", "10", ",", "24"}], "}"}]}], ";"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.509597551484375*^9, 3.509597558140625*^9}, 
   3.509598835375*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"ddiff", "=", 
   RowBox[{
    RowBox[{"Reap", "[", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", 
       RowBox[{"i", "<=", 
        RowBox[{
         RowBox[{"Length", "[", "ptest", "]"}], "-", "1"}]}], ",", 
       RowBox[{"i", "++"}], ",", "\n", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"j", "=", 
          RowBox[{"i", "+", "1"}]}], ",", 
         RowBox[{"j", "<=", 
          RowBox[{"Length", "[", "ptest", "]"}]}], ",", 
         RowBox[{"j", "++"}], ",", "\n", 
         RowBox[{
          RowBox[{"Sow", "[", 
           RowBox[{"EuclideanDistance", "[", 
            RowBox[{
             RowBox[{"ptest", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"ptest", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], "]"}], "]"}], ";"}]}], "\n", 
        "]"}]}], "]"}], "\n", "]"}], "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ";"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.50959866803125*^9, 3.50959876234375*^9}, 
   3.509598835375*^9}],

Cell[BoxData[
 RowBox[{"Min", "[", "ddiff", "]"}]], "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.509598765265625*^9, 3.50959878178125*^9}, 
   3.509598835375*^9}],

Cell[BoxData["1.5211491332608467`"], "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.},
 CellChangeTimes->{{3.509598767046875*^9, 3.509598782296875*^9}, 
   3.509598835375*^9}]
}, Closed]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"geteer", "[", 
  RowBox[{
   RowBox[{"Partition", "[", 
    RowBox[{"stpc11", ",", "2"}], "]"}], ",", "tg", ",", "classifyDist1"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.50890544359375*^9, 3.508905445640625*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"5.354999999999908`", ",", 
   RowBox[{"{", 
    RowBox[{"0.1`", ",", "0.1`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.508905478390625*^9}]
}, Open  ]]
},
AutoGeneratedPackage->Automatic,
WindowSize->{1016, 520},
WindowMargins->{{0, Automatic}, {0, Automatic}},
ShowSelection->True,
FrontEndVersion->"8.0 for Microsoft Windows (32-bit) (November 7, 2010)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 146, 3, 39, "Input"],
Cell[716, 27, 594, 12, 25, "Message"],
Cell[1313, 41, 593, 12, 28, "Message"],
Cell[1909, 55, 593, 12, 25, "Message"],
Cell[2505, 69, 682, 13, 22, "Message"]
}, Open  ]],
Cell[3202, 85, 581, 12, 94, "Input"],
Cell[CellGroupData[{
Cell[3808, 101, 409, 10, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[4220, 113, 2039, 37, 148, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 InitializationCell->True],
Cell[6262, 152, 7868, 192, 715, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 InitializationCell->True],
Cell[14133, 346, 1865, 30, 94, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.},
 InitializationCell->True]
}, Open  ]],
Cell[CellGroupData[{
Cell[16035, 381, 474, 10, 67, "Input"],
Cell[16512, 393, 364, 6, 38, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16913, 404, 630, 15, 67, "Input"],
Cell[17546, 421, 350, 7, 38, "Output"]
}, Open  ]],
Cell[17911, 431, 574, 16, 94, "Input"],
Cell[18488, 449, 198, 5, 39, "Input"],
Cell[CellGroupData[{
Cell[18711, 458, 152, 3, 39, "Input"],
Cell[18866, 463, 96, 1, 38, "Output"]
}, Closed]],
Cell[CellGroupData[{
Cell[18999, 469, 221, 5, 31, "Input"],
Cell[19223, 476, 3001, 40, 289, "Output"]
}, Closed]],
Cell[CellGroupData[{
Cell[22261, 521, 343, 7, 59, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[22607, 530, 2020, 43, 148, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}]
}, Closed]],
Cell[CellGroupData[{
Cell[24664, 578, 746, 13, 59, "Input"],
Cell[25413, 593, 206, 4, 38, "Output"]
}, Open  ]],
Cell[25634, 600, 867, 24, 67, "Input"],
Cell[26504, 626, 174, 3, 39, "Input"],
Cell[CellGroupData[{
Cell[26703, 633, 279, 7, 39, "Input"],
Cell[26985, 642, 163, 2, 38, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27185, 649, 198, 5, 39, "Input"],
Cell[27386, 656, 7016, 203, 364, "Output"]
}, Open  ]],
Cell[34417, 862, 336, 9, 39, "Input"],
Cell[34756, 873, 296, 9, 39, "Input"],
Cell[35055, 884, 247, 5, 39, "Input"],
Cell[35305, 891, 92, 1, 39, "Input"],
Cell[35400, 894, 745, 19, 121, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[36148, 915, 116, 3, 38, "Output"],
Cell[CellGroupData[{
Cell[36289, 922, 570, 16, 67, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[36862, 940, 763, 16, 256, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[37628, 958, 254, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[37885, 964, 252, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[38140, 970, 255, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[38398, 976, 256, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[38657, 982, 260, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[38920, 988, 255, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[39178, 994, 255, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[39436, 1000, 256, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[39695, 1006, 1399, 44, 202, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}]
}, {11}]],
Cell[CellGroupData[{
Cell[41129, 1055, 1592, 47, 283, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[42724, 1104, 478, 9, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[43205, 1115, 468, 10, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}]
}, Closed]],
Cell[CellGroupData[{
Cell[43710, 1130, 650, 13, 59, "Input"],
Cell[44363, 1145, 365, 8, 38, "Output"]
}, Closed]],
Cell[44743, 1156, 364, 10, 31, "Input"],
Cell[CellGroupData[{
Cell[45132, 1170, 119, 2, 39, "Input"],
Cell[45254, 1174, 119, 2, 38, "Output"]
}, Closed]],
Cell[CellGroupData[{
Cell[45410, 1181, 1053, 27, 59, "Input"],
Cell[46466, 1210, 540, 12, 38, "Output"]
}, Closed]],
Cell[47021, 1225, 2342, 57, 248, "Input"],
Cell[CellGroupData[{
Cell[49388, 1286, 579, 15, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[49970, 1303, 2516, 46, 435, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}]
}, Closed]],
Cell[CellGroupData[{
Cell[52523, 1354, 330, 7, 59, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[52856, 1363, 500, 14, 94, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[53359, 1379, 227, 5, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[53589, 1386, 543, 15, 67, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}]
}, Closed]],
Cell[CellGroupData[{
Cell[54169, 1406, 1329, 43, 194, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10002.}],
Cell[55501, 1451, 556, 16, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10002.}],
Cell[56060, 1469, 2395, 45, 408, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10002.}]
}, Closed]],
Cell[CellGroupData[{
Cell[58492, 1519, 1947, 54, 302, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10003.}],
Cell[60442, 1575, 487, 9, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10003.}],
Cell[60932, 1586, 224, 5, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10003.}],
Cell[61159, 1593, 376, 7, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10003.}]
}, Closed]],
Cell[61550, 1603, 2996, 74, 518, "Input"],
Cell[CellGroupData[{
Cell[64571, 1681, 806, 18, 94, "Input"],
Cell[65380, 1701, 240, 5, 38, "Output"]
}, Open  ]],
Cell[65635, 1709, 1594, 38, 94, "Input"],
Cell[67232, 1749, 231, 5, 39, "Input"],
Cell[CellGroupData[{
Cell[67488, 1758, 234, 4, 67, "Input"],
Cell[67725, 1764, 136, 3, 38, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[67898, 1772, 498, 11, 39, "Input"],
Cell[68399, 1785, 2053, 38, 404, "Output"]
}, Closed]],
Cell[CellGroupData[{
Cell[70489, 1828, 335, 7, 31, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10004.}],
Cell[70827, 1837, 367, 9, 67, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10004.}],
Cell[71197, 1848, 248, 5, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10004.}]
}, Closed]],
Cell[CellGroupData[{
Cell[71482, 1858, 1674, 44, 194, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[73159, 1904, 188, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[73350, 1910, 187, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[73540, 1916, 534, 14, 67, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[74077, 1932, 191, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[74271, 1938, 1501, 41, 202, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}]
}, Closed]],
Cell[75787, 1982, 808, 24, 59, "Input"],
Cell[CellGroupData[{
Cell[76620, 2010, 428, 10, 39, "Input"],
Cell[77051, 2022, 3147, 55, 374, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[80235, 2082, 585, 14, 94, "Input"],
Cell[80823, 2098, 1304, 20, 139, "Output"],
Cell[82130, 2120, 138, 2, 38, "Output"]
}, Open  ]],
Cell[82283, 2125, 264, 8, 67, "Input"],
Cell[CellGroupData[{
Cell[82572, 2137, 368, 7, 39, "Input"],
Cell[82943, 2146, 111, 1, 38, "Output"]
}, Open  ]],
Cell[83069, 2150, 213, 5, 39, "Input"],
Cell[83285, 2157, 1523, 39, 148, "Input"],
Cell[CellGroupData[{
Cell[84833, 2200, 183, 4, 39, "Input"],
Cell[CellGroupData[{
Cell[85041, 2208, 66, 1, 29, "Print"],
Cell[85110, 2211, 67, 1, 29, "Print"],
Cell[85180, 2214, 68, 1, 29, "Print"],
Cell[85251, 2217, 69, 1, 29, "Print"],
Cell[85323, 2220, 69, 1, 29, "Print"],
Cell[85395, 2223, 68, 1, 29, "Print"],
Cell[85466, 2226, 69, 1, 29, "Print"],
Cell[85538, 2229, 68, 1, 29, "Print"],
Cell[85609, 2232, 67, 1, 29, "Print"],
Cell[85679, 2235, 70, 1, 29, "Print"],
Cell[85752, 2238, 70, 1, 29, "Print"],
Cell[85825, 2241, 70, 1, 29, "Print"],
Cell[85898, 2244, 67, 1, 29, "Print"],
Cell[85968, 2247, 69, 1, 29, "Print"],
Cell[86040, 2250, 70, 1, 29, "Print"],
Cell[86113, 2253, 69, 1, 29, "Print"],
Cell[86185, 2256, 69, 1, 29, "Print"],
Cell[86257, 2259, 70, 1, 29, "Print"],
Cell[86330, 2262, 70, 1, 29, "Print"],
Cell[86403, 2265, 69, 1, 29, "Print"],
Cell[86475, 2268, 67, 1, 29, "Print"],
Cell[86545, 2271, 70, 1, 29, "Print"],
Cell[86618, 2274, 70, 1, 29, "Print"],
Cell[86691, 2277, 70, 1, 29, "Print"],
Cell[86764, 2280, 70, 1, 29, "Print"],
Cell[86837, 2283, 69, 1, 29, "Print"],
Cell[86909, 2286, 69, 1, 29, "Print"],
Cell[86981, 2289, 70, 1, 29, "Print"],
Cell[87054, 2292, 70, 1, 29, "Print"],
Cell[87127, 2295, 70, 1, 29, "Print"],
Cell[87200, 2298, 70, 1, 29, "Print"],
Cell[87273, 2301, 69, 1, 29, "Print"],
Cell[87345, 2304, 69, 1, 29, "Print"],
Cell[87417, 2307, 70, 1, 29, "Print"],
Cell[87490, 2310, 70, 1, 29, "Print"],
Cell[87563, 2313, 67, 1, 29, "Print"],
Cell[87633, 2316, 70, 1, 29, "Print"],
Cell[87706, 2319, 68, 1, 29, "Print"],
Cell[87777, 2322, 69, 1, 29, "Print"],
Cell[87849, 2325, 70, 1, 29, "Print"],
Cell[87922, 2328, 69, 1, 29, "Print"],
Cell[87994, 2331, 70, 1, 29, "Print"],
Cell[88067, 2334, 65, 1, 29, "Print"],
Cell[88135, 2337, 68, 1, 29, "Print"],
Cell[88206, 2340, 70, 1, 29, "Print"],
Cell[88279, 2343, 70, 1, 29, "Print"],
Cell[88352, 2346, 70, 1, 29, "Print"],
Cell[88425, 2349, 70, 1, 29, "Print"],
Cell[88498, 2352, 67, 1, 29, "Print"],
Cell[88568, 2355, 69, 1, 29, "Print"],
Cell[88640, 2358, 70, 1, 29, "Print"],
Cell[88713, 2361, 70, 1, 29, "Print"],
Cell[88786, 2364, 70, 1, 29, "Print"],
Cell[88859, 2367, 70, 1, 29, "Print"],
Cell[88932, 2370, 69, 1, 29, "Print"],
Cell[89004, 2373, 64, 1, 29, "Print"],
Cell[89071, 2376, 68, 1, 29, "Print"],
Cell[89142, 2379, 70, 1, 29, "Print"],
Cell[89215, 2382, 69, 1, 29, "Print"],
Cell[89287, 2385, 70, 1, 29, "Print"],
Cell[89360, 2388, 68, 1, 29, "Print"],
Cell[89431, 2391, 69, 1, 29, "Print"],
Cell[89503, 2394, 70, 1, 29, "Print"],
Cell[89576, 2397, 67, 1, 29, "Print"],
Cell[89646, 2400, 67, 1, 29, "Print"],
Cell[89716, 2403, 68, 1, 29, "Print"],
Cell[89787, 2406, 70, 1, 29, "Print"],
Cell[89860, 2409, 70, 1, 29, "Print"],
Cell[89933, 2412, 70, 1, 29, "Print"],
Cell[90006, 2415, 64, 1, 29, "Print"],
Cell[90073, 2418, 70, 1, 29, "Print"],
Cell[90146, 2421, 70, 1, 29, "Print"],
Cell[90219, 2424, 69, 1, 29, "Print"],
Cell[90291, 2427, 64, 1, 29, "Print"],
Cell[90358, 2430, 70, 1, 29, "Print"],
Cell[90431, 2433, 68, 1, 29, "Print"],
Cell[90502, 2436, 68, 1, 29, "Print"],
Cell[90573, 2439, 70, 1, 29, "Print"],
Cell[90646, 2442, 69, 1, 29, "Print"],
Cell[90718, 2445, 70, 1, 29, "Print"],
Cell[90791, 2448, 70, 1, 29, "Print"],
Cell[90864, 2451, 69, 1, 29, "Print"],
Cell[90936, 2454, 70, 1, 29, "Print"],
Cell[91009, 2457, 70, 1, 29, "Print"],
Cell[91082, 2460, 68, 1, 29, "Print"],
Cell[91153, 2463, 70, 1, 29, "Print"],
Cell[91226, 2466, 69, 1, 29, "Print"],
Cell[91298, 2469, 67, 1, 29, "Print"],
Cell[91368, 2472, 70, 1, 29, "Print"],
Cell[91441, 2475, 69, 1, 29, "Print"],
Cell[91513, 2478, 69, 1, 29, "Print"],
Cell[91585, 2481, 70, 1, 29, "Print"],
Cell[91658, 2484, 70, 1, 29, "Print"],
Cell[91731, 2487, 70, 1, 29, "Print"],
Cell[91804, 2490, 68, 1, 29, "Print"],
Cell[91875, 2493, 70, 1, 29, "Print"],
Cell[91948, 2496, 70, 1, 29, "Print"],
Cell[92021, 2499, 70, 1, 29, "Print"],
Cell[92094, 2502, 70, 1, 29, "Print"]
}, Closed]],
Cell[92179, 2506, 98, 1, 36, "Output"]
}, Open  ]],
Cell[92292, 2510, 243, 7, 39, "Input"],
Cell[CellGroupData[{
Cell[92560, 2521, 378, 8, 67, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[92941, 2531, 728, 16, 67, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[93672, 2549, 2022, 37, 464, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}]
}, Closed]],
Cell[95709, 2589, 991, 21, 86, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[CellGroupData[{
Cell[96725, 2614, 960, 29, 94, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[97688, 2645, 248, 5, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[97939, 2652, 336, 7, 67, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[98278, 2661, 240, 5, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[98521, 2668, 249, 6, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[98773, 2676, 466, 10, 67, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[99242, 2688, 206, 4, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[99451, 2694, 191, 4, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[99645, 2700, 290, 7, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[99938, 2709, 285, 5, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[100226, 2716, 496, 13, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[100725, 2731, 783, 23, 121, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[101511, 2756, 368, 9, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[101882, 2767, 1143, 32, 148, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[103028, 2801, 201, 4, 39, "Input",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}],
Cell[103232, 2807, 188, 3, 38, "Output",
 CellGroupingRules->{GroupTogetherGrouping, 10001.}]
}, Closed]],
Cell[CellGroupData[{
Cell[103457, 2815, 240, 6, 31, "Input"],
Cell[103700, 2823, 187, 5, 38, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
